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In tropospheric scatter transmission beyond the horizon, the amplitude, 
phase and frequency of a received sine wave exhibit random fluctuations 
owing to variable multipath transmission and noise. The probability of 
errors in digital transmission over such random multipath media has been 
dealt with in the literature on the premise of flat Rayleigh fading over the 
band occupied by the spectrum of transmitted pidses. This is a legitimate 
approximation at low trans?nission rates, such that the pulse spectrum is 
adequately narrow, bid not at high digital transmission rates. The proba- 
bility of errors is determined here also for high transmission rates, such that 
selective fading over the pulse spectrum band must be considered. Such 
selective fading gives rise to pxdse distortion and resultant intersymbol 
interference that may cause errors even in the absence of noise. 

Troposcatter transmission can be approximated by an idealized multi- 
path model in which the amplitudes of signal wave components received 
over different paths vary at random and in which there is a linear variation 
in transmission delay with a maximum departure ±A from the mean 
delay. Various statistical transmission parameters are determined on this 
premise, among them the probability distnbidion of amplitude and phase 
fluctuations and of derivatives thereof with respect to time and with respect 
to frequency. The probability of errors in the absence of noise owing to such 
fluctuations is determined together with the probability of errors owing to 
noise, for digital transmissio?i by binary PM and FM. Charts are pre- 
sented, from which can be determined the combined probability of errors from 
various sources, as related to the ti'a?ismission rate and certain basic param- 
eters of troposcatter links. 
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INTRODUCTION 

In tropospheric transmission beyond the horizon, narrow-beam 
transmitting and receiving antennas are used in a frequency range from 
about 400 to 10,000 megacycles. The received wave can be considered 
the sum of a large number of components of varying amplitudes, re- 
sulting from a multiplicity of reflections within the common volume at 
the intersection of the antenna beams. These various components arrive 
with different transmission delays owing to path-length differences, 
and each will exhibit a variation in amplitude owing to structural 
changes within the common volume, caused largely by winds. When a 
steady -state sine wave is transmitted, the received wave will conse- 
quently exhibit variations in its envelope and phase, commonly referred 
to as fading. When a signal wave is transmitted, its various frequency 
components will suffer unwanted amplitude and phase variations with 
resultant transmission impairments that depend on the particular 
carrier modulation method. These impairments are discussed herein 
for digital transmission by carrier phase and frequency modulation. 

Various properties of the transmittance of troposcatter channels 
have been dealt with in several publications. 1 2 ,3 - 4 These properties 
include the expected average path loss and systematic seasonal varia- 
tions from the average, together with the probability distributions of 
slow and rapid fading or fluctuations from the mean. Other important 
properties from the standpoint of systems design and performance are 
the distribution of duration of fades and the fading rapidity or rate. 

The above various properties relate to transmittance variations with 
time at a particular frequency. Of basic importance is also the variation 
in transmittance with frequency at any instant, i.e., the amplitude and 
phase characteristics of trophospheric channels. These will be highly 
variable quantities, as illustrated in Fig. 1. At a fixed instant the 
characteristics may be as indicated in Fig. 1(a) and at a later instant 
as in Fig. 1(b). Such fluctuations will give rise to a distortion of the 
spectrum of received signals, with resultant transmission impairments 
of various kinds, depending on the modulation method. In addition, 
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random noise at the receiver input must be considered as in conventional 
stable channels. Owing to the above random fluctuations, diversity 
transmission is ordinarily required to insure adequate performance. 

At present, frequency modulation is used for transmission of multi- 
plexed voice channels over troposcatter links. With this method, pro- 
nounced intermodulation noise is encountered 5 ' 6 owing to the types of 
transmittance variations with frequency indicated in Fig. 1. With digi- 
tal transmission, these variations will give rise to pulse distortion and 
resultant intersymbol interference that may severely limit the trans- 
mission rate. 

In evaluation of error probabilities in digital transmission, it is neces- 
sary to consider variations in the average path loss over a convenient 
period, such as an hour, relative to the average over a much longer 
period, say a month. These slow fluctuations in loss are closely approxi- 
mated by the log-normal law, i.e., the loss in db follows the normal 
law. 1 In addition, consideration must be given to rapid fluctuations in 
loss relative to the above hourly averages. These are closely approxi- 
mated by the Rayleigh law, which also applies for the envelope of 
narrow-band random noise. They are ordinarily more important than 
slow fluctuations, particularly in digital transmission, in that they cannot 
be fully compensated for by automatic gain control. Nearly all theoreti- 
cal analyses of error probabilities in digital transmission over fading 
channels are based on a Rayleigh distribution together with various 
other simplifying assumptions, as outlined below. 

The simplest assumption is flat or nonselective Rayleigh fading over 
the channel band, in conjunction with a sufficiently slow fading rate 
such that changes over a few pulse intervals can be disregarded. These 
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Fig. 1 — Illustrative variations in attenuation and phase characteristics with 
frequency at two instants ti and tt . 
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are legitimate premises in transmission over line-of-sight radio links, 
where fading is much slower than on tropospheric links and is virtually 
nonselective over rather wide bands. With these simplifying assumptions 
Turin 7 has determined error probabilities in binary transmission over 
noisy channels with ideal synchronous (coherent) detection and envelope 
(noncoherent) detection. His analysis includes the effect of correlation 
between successive pulses and also postulates a nonfading signal com- 
ponent, such that the results in one limit also apply for nonfading chan- 
nels. 

On the same premise of slow, flat Rayleigh fading, Pierce 8 has deter- 
mined the optimum theoretical diversity improvement for frequency 
shift keying with dual filter reception employing coherent and non- 
coherent detection of the filter outputs. Dual filter detection is ordi- 
narily assumed in place of the usual method of frequency discriminator 
detection that does not lend itself as readily to theoretical analysis. 

The error probability with two-phase and four-phase modulation 
with differential phase detection has been determined by Voelcker 9 on 
the premise of flat Rayleigh fading at such a rate that the change in 
phase over a pulse interval must be considered. Moreover, he considers 
the probability of both single and double digital errors, with both single 
and dual diversity transmission. 

Voelcker's analysis is applicable to transmission at a sufficiently slow 
rate such that amplitude and phase distortion can be ignored over the 
relatively narrow band of the pulse spectra. However, it does not apply 
to high-speed digital transmission that requires sufficiently wide pulse 
spectra such that the amplitude and phase distortion indicated in Fig. 
1 must be considered. For this case the duration of pulses will be so 
short that the phase changes considered by Voelcker can be disregarded. 
Instead, it now becomes necessary to take into account pulse distortion 
and resultant intersymbol interference caused by the erratic variations 
with frequency in the amplitude and phase characteristics illustrated 
in Fig. 1. An evaluation is made herein of error probabilities on the 
latter account, which has not been considered in previous publications.* 

From the solutions for the above two limiting cases of low and high 
transmission rates, it is possible by simple graphical methods to esti- 
mate the error probability for the general case in which both time and 
frequency variations in the amplitude and phase characteristics must 
be considered. Charts are presented of error probabilities in digital 
transmission by binary PM and FM as related to various basic param- 
eters of tropospheric scatter links and of the signals. Among these 



* For reference to a recent related paper, see Section 8.9. 
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parameters are the average signal-to-noise ratio, the bandwidth of the 
pulse spectrum, the fading bandwidth of the troposcatter link, and the 
maximum departure from the mean transmission delay, which is re- 
lated to the length of the link and the antenna beam angles. 

The analysis shows that a principal source of pulse distortion and 
resultant transmission impairments is a component of quadratic phase 
distortion. On this premise, an evaluation has been made in a companion 
paper* of intermodulat ion distortion in analog transmission by FM and 
PM, that conforms well with the results of measurements. 5 - 6 

I. CHANNEL TRANSMISSION CHARACTERISTICS 

1.1 General 

Transmission performance with any modulation method depends on 
the statistical properties of the signals and of channel noise, together 
with various properties of the channel transmittance or transmission- 
frequency characteristic. When the latter varies with time, the usual 
methods of determining network response to specified input waves must 
be modified in various respects, that result in appreciable complications 
in the analytical methods 10 and in certain conceptual difficulties. How- 
ever, when the time variations in transmittance are slow in relation to 
those in the input waves, it is legitimate to assume that the trans- 
mittances are constant over an appreciable number of pulse intervals. 
With relatively slow random fluctuations as encountered in troposcatter 
systems at representative transmission rates, it is thus permissible to 
determine the responses for various essentially time invariant transmit- 
tances that can be encountered. In evaluating transmission performance, 
the various transmittances that can be encountered must be weighted 
or averaged statistically in a manner that depends on the signal prop- 
erties and the modulation method. 

Among the statistical properties of troposcatter transmittances are 
the probability distribution of the envelope of received carrier waves 
together with the autocorrelation function of the envelope with re- 
spect to time and with respect to frequency. These are discussed here, 
while other statistical properties will be considered in later sections. 

1.2 Tropospheric Scatter Waves 

To determine an appropriate model for the random process in trop- 
ospheric scatter transmission, it is necessary to consider the physics 

* See part 2 of this issue of the B.S.T.J., to appear. 
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of this phenomenon, as dealt with in various publications. Though these 
may differ in their assumptions regarding the exact mechanism of the 
reflections, they appear to agree that they occur as a result of hetero- 
geneities within the common antenna volume indicated in Fig. 2. If the 
transmission medium were uniform, no reception would be possible. 
Owing to the numerous heterogeneities in the common volume, a very 
large number of reflections will occur, and the received wave can be 
considered the sum of a large number of components of different ampli- 
tudes and different transmission delays. Over any short interval, the 
envelope of a received sine wave will depend on the frequency, as will 
the phase. Because of variations in the heterogeneities caused largely 
by winds, the envelope and phase of a received carrier will vary with 
time. 



Fig. 2 — Illustrative antenna beams and common antenna volume. 

The transmittance of troposcatter channels is dealt with here, based 
on an idealized model discussed further in the Appendix, and certain 
statistical parameters obtained from experimental data are discussed. 
Two limiting cases that permit simplified analysis are considered. In 
one case the transmission band is assumed sufficiently narrow, such that 
the attenuation characteristic can be considered constant and the phase 
characteristic linear over the narrow band. There will then be fluctua- 
tions with time in the attenuation accompanied by independent varia- 
tions in the slope of the phase characteristic, a condition referred to 
as nonselective flat fading and ordinarily assumed in random multipath 
digital transmission theory. The other limiting case is that of digital 
transmission at a sufficiently high rate so that time variations in the 
transmittance can be disregarded over an appreciable number of pulse 
intervals. In this case it is necessary to consider erratic variations with 
frequency in both the attenuation and phase characteristics. 
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Fig. 3 — Illustrative dependence of envelope and phase of transmittance with 
frequency u from a reference frequency wo at a specified time t\ . 

1.3 Troposcatter Transmittance 

Let a sine wave of frequency oj be transmitted, and let w = a- + u, 
as indicated in Fig. 3, where o> is a conveniently chosen reference fre- 
quency. In complex notation the received wave is then of the general 
form 

e(u,t) = r(u,t) exp[— i<p(u,t)] exp(uat) (1) 

where r(u,t) and <p(u,t) are random variables of the time t for a fixed 
to or u, and of u for a fixed time t. The channel transmittance is then 

T(u,t) = r(u,t) exp[-i<p(u,t)]. (2) 

The following general relations apply 

r(u,t) = [U\u,t) + V\u,t))> 

<p(u,l) = tan -1 [V(u,t)/U(u,t)]. 



(3) 
(4) 



As shown in the Appendix, in the case of idealized tropospheric 
channels the functions U and V can be represented in the following form 



U(ujt) = £ aj(t] 



sin (jt — u A) 
Jt — uA 



(5) 
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V(u,l) = ± 6,(0 Si "° V ~ A UA) (6) 

/=_« Jtt — UA 



where 



A = maximum departure from mean transmission delay 
owing to path length differences. 

In (5) and (G) the coefficients a.j(t) and bj(l) vary at random with 
time t and for a given t vary at random with j. Owing to the latter 
variation with j, there will be a random variation in U and V with the 
frequency u taken in relation to the reference frequency wo . 

Equations for an idealized troposcatter channel, as given in the 
Appendix, show that oy(0 is related to the sum A(x,t) + A(—x,t) of 
two random processes and bj(t) to the difference A(x,t) — A(—x,t). 
The two random processes A(x,t) and A{— x,t) will have equal rms 
amplitudes, in which case a,-(<) and bj(t) will have zero correlation 
coefficient. They will then also be independent random variables, pro- 
vided A(x,t) and A(—x,t) have a Gaussian probability distribution, 
which appeal's to be a legitimate approximation since each will be the 
sum of waves from a large number of reflections. 

A further assumption underlying (5) and (G) is that there is an in- 
finite number of transmission paths. An additional approximation that 
will be made in the following analysis is that there will be independent 
random fluctuations in the signal components received over the various 
paths. Actually there will be some correlation between the fluctuations, 
particularly for paths with small separation. In effect, there will be a 
limited number of essentially independently fading paths. 

The above assumptions entail certain statistical properties of tropo- 
scatter channels, as outlined below for time and frequency variations. 

1.4 Transmission Loss Fluctuations 

On troposcatter links there is a certain average transmission loss over 
a year, which depends on the length of the link, on the properties of the 
terrain and on climatic conditions. Experimental data indicate that 
there will be systematic monthly and seasonal departures from this 
yearly average, owing principally to slow temperature changes. The 
average loss during a winter month may thus be up to 20 db greater 
than the average during a summer month. That is, the departure in 
transmission loss from the yearly mean may be ±10 db. 

During each month there will be a more or less random fluctuation 
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in the hourly average loss from the mean of the month. This fluctuation 
has been found to be almost independent of frequency and seems to be 
associated with the variations in average refraction of the atmosphere 
and resultant variation in the bending of beams. This fluctuation in 
the hourly average loss relative to the monthly average has been found 
to follow closely the log-normal law. That is to say, let the monthly 
median loss be 

a m = —In f OT 2 (7) 

and the hourly average loss be 

a = —In f (8) 

where In = log £ , f m is the monthly rms amplitude of the envelope 
r(u,t), and f the rms amplitude over an hour. (Other reference times 
could have been chosen, as will appear below.) 

The probability that the average hourly loss exceeds a specified value 
«i = In fi is then given by 



P(« ^ «0 = fi 1 - erf -^^ J 



(9) 



where erf is the error function and <r a the standard deviation in trans- 
mission loss expressed in nepers, when a and a m are expressed in nepers 
as above. For links 100 to 200 miles in length, a representative value of 
a a appears to be about 0.9 neper (8 db). 

In addition to the above slow variations in the average hourly loss, 
there will be more rapid fluctuations in the envelope r(u,t), owing to 
changes in the multipath transmission structure caused principally by 
winds. This type of fluctuation follows a Rayleigh distribution law. 
According to this law the probability that the instantaneous value r 
of the envelope exceeds a specified value n is 

P(r> n) = exp(-r 1 2 /f 2 ) (10) 

where f is the hourly rms value referred to above. 

It may be noted that while the log-normal law for slow variation has 
been determined solely by measurements, the Rayleigh law for rapid 
fluctuations follows by theory when the received wave is the sum of a 
large number of variable components. 

The probability distribution (10) can be related to the monthly rms 
value of r(u,t) with the aid of (9) by 

P(r > r,) = f p(f) exp(-r 1 2 /r 2 ) df (11) 
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where p(f) is the probability density function corresponding to (9), 
which is 

P(r) = ^7=-: exp [ - [In f 2 /rjf/2* a 2 }. ( 12) 

It will be recognized that (11) will yield the same result regardless 
of the period over which the rms value f is taken, since r simply plays 
the role of an intermediate parameter that disappears after integration. 

The above probability functions relating to average loss or the dis- 
tribution of the instantaneous values of r(u,t), are independent of the 
frequency. In addition to the above distribution there are others which 
are important from the standpoint of transmission systems design and 
performance, as discussed in the following section. 

1.5 Time Autocorrelation Functions of Trans mil tance 

Expressions for the probabilities of rapid changes in the amplitude 
and phase of the transmittance with time will be considered in Section 
II. These involve the autocorrelation functions of the components U 
and V defined by (5) and (6), or the corresponding power spectra. Both 
have the same autocorrelation function and power spectrum, so that 
only U(u,t) needs to be considered. 

The time autocorrelation function of U{u,t) depends on the variation 
in aj(t) with time. These are related to changes in the physical structure 
of the common volume and to resultant variations in the heterogeneities 
that are responsible for tropospheric transmission. The rate at which 
these occur depends on the velocity and directions of winds and on 
temperature changes. Under these conditions the autocorrelation func- 
tion will vary with time, and it becomes necessary to consider a certain 
median autocorrelation function and corresponding power spectrum, 
as discussed in Section 1.6. 

Let ^(r) be the autocorrelation function of variations in U(u,t) with 
/. The corresponding one-sided power spectrum is then 

9 r 

W(y) = - / tf(r) cosyrdr (13) 

7T •'0 

where y is used to designate the radian frequency of spectral compo- 
nents to avoid confusion with the frequency co of the transmitted wave. 
The autocorrelation function \Ht) or the corresponding power spec- 
trum W(y) of the components £7 and V cannot be determined as readily 
by measurements as the autocorrelation function & t (t) of the envelope. 
The latter is related to ^(t) by" 
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* r ( T ) = *(0){2E[k(t)] - [1 - k\t)]K[k(t)]} (14) 

where 

k(t) = ¥(t)/¥(0) (15) 

E = complete elliptic integral of second kind 
K = complete elliptic integral of first kind. 

For t = 0, ^r(O) = 2^(0). Hence the autocorrelation coefficient of 
the envelope can be written 

Kr (r) - E[k(t)] - i[l - k 2 (t)]K[k(t)]. (16) 

With the aid of (16), the autocorrelation coefficient k(t) of each 
quadrature component can be determined from measurements of <c r (r). 

1.6 Observed Time Autocorrelation 

Observations of the autocorrelation function of rapid fluctuations 
indicate that the autocorrelation function ¥(t) of the components U 
and V is nearly Gaussian and is given by 

*( T ) = *(0)exp(-<rV/2). (17) 

The corresponding power spectrum obtained from (13) is 

W(y) = *(0)(2/tt(7 2 ) J exp(-y 2 /2a 2 ) (18) 

where ^(0) is the average power in each component as obtained with 
t = Oin (17) . 

The equivalent bandwidth of a fiat power spectrum W(y) = W(0) 
is given by 

y = VUM^ — 1.25*. (19) 

As noted in Section 1.5, there will be a certain median autocorrelation 
function and corresponding median values of the power spectrum, of 
a and of y. Measurements 2 indicate that these median values depend on 
the antenna beamwidths and that the fading rate is not quite propor- 
tional to frequency. Furthermore, there can be appreciable departure 
from the median values. From measurements of the median number of 
fades per minute, the median value of a can be determined, with the aid 
of equation (26) in Ref. 2. These measurements indicate that for a 
particular antenna arrangement a & 0.1 cps at 460 mc and about 1.3 
cps at 4110 mc. The corresponding equivalent bandwidths of a flat 
power spectrum are thus y 9H 0.125 cps, or 0.8 radian/sec. at 460 mc, 
and a « 1.6 cps, or about 10 radians/sec. at 4110 mc. The measurements 
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further indicate that there is a probability of about 0.01 that the fading 
rate exceeds the median value by a factor of about 7 at 460 mc and a 
factor of about 3.5 at 4110 mc. 

1.7 Frequency Correlation Function of Transmittance 

Returning to (5) and (G), let the time t be fixed, and consider varia- 
tions in U and V with u. The coefficients a, and b } will then have certain 
values that vary with j, and there will be a certain variation in U and 
V with u. At a different time there will be another set of coefficients 
and a different variation with u. The form of (5) and (0) indicates that 
if u is regarded as a time variable and A as a frequency, U(u) would 
be the variation in time owing to impulses of amplitudes a, and bj 
impinging at time intervals ir on a flat low-pass filter of bandwidth A. 
That is to say, the autocorrelation function of components U and V 
for a difference v = coo — wi in frequency is 

*(r) = ¥(0)(sin "A/j/A). (20) 

The corresponding power spectrum of the variation in U and V with 
frequency 8 is 

2 f°° 

W(8) =- / *0) cos v8 dv (21) 

IT Jf) 

= *(0) for < 5 < A 

(22) 
= for A < 8. 

When V(v) is given, it is possible to determine the autocorrelation 
function ^ r (v) for variations in r(u,t) with u. Expression (14) applies 
with v in place of r, for the autocorrelation function of time variation 
with frequency. 

For an autocorrelation function (20), the corresponding correlation 
coefficient is 

k(v) = (sin fA/VA). (23) 

The corresponding autocorrelation coefficient of the envelope, as ob- 
tained from (10), is 

/ v 7-,/sin vA\ if, sin 2 vA~] ,Vsin uA\ ,_.v 

For various values of vA the correlation function of the envelope is 
given in Table I and is shown in Fig. 4. 
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Table I — Autocorrelation Function of Envelope 
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The autocorrelation functions (23) and (24) apply for certain idealized 
conditions outlined in the Appendix and in Section 1.3. For one thing, 
the average power received over each elementary path is assumed the 
same. For another, a linear variation in the transmission delay with 
angular deviation from the mean paths is assumed, with maximum 
departures ±A from the mean delay. Furthermore, an infinity of trans- 
mission paths is assumed, with independent random fluctuations in the 
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Fig. 4 — Frequency autocorrelation coefficient k,(») of envelope for autocorre- 
lation coefficient k(») of components U and V. 

signal components received over the various paths, though there will 
be some correlation between the fluctuations in the signal components 
received over various paths. 

In spite of the various approximations, it appears possible to obtain 
a reasonably satisfactory conformance with the results of measure- 
ments of the autocorrelation functions of the envelope, as shown in 
Section 1.9. 



DIGITAL TROPOSCATTER 



157 



1.8 Differential Transmission Delay -\ 

Exact determination of the equivalent maximum departure from the 
mean transmission delay requires consideration of the beam patterns 
as affected by scattering. On the approximate basis of equivalent beam 
angles a, the following relation applies, with notation as indicated in 
Fig. 5 



La + B 
v 2 



+ 



a + B 

■2 



(25) 



where 3 ^ a, v is the velocity of propagation in free space, L is the 
length of the link, and 



L 

2R 



L 



2RoK 



(26) 



where 7? is the radius of the earth and the factor K is ordinarily taken 
as 4/3. 

The equivalent beam angle a from midbeam to the 3-db loss point 
depends on the free-space antenna beam angle a and on the effect of 
scatter, which is related in a complex manner to ao and the length L, 
or alternately 6. Narrow-beam antennas as now used in actual systems 
are loosely defined by a ^ 26/3. For these a £& <x on shorter links, 
while on longer links a > «o owing to beam-broadening by scatter. 
Analytical determination of a for longer links appears difficult, and only 




Fig. 5 — Definition of antenna beam angles a, take-off angle /3 and chord angle 
6 to midbeam. With different angles at the two ends, the mean angles are used in 
expressions for A. In applications to actual beams, a would be the angle to the 
3-db loss point. 
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limited experimental data are available at present. For broad-beam 
antennas, a > 20/3 and beam-broadening by scatter is in theory in- 
appreciable. 

By way of numerical example, let L = 170 miles and K = 4/3, in 
which case = 0.016 radian. Since a = 0.004 radian « 20/3, it is 
permissible to take a = a Q . With « <x , (25) gives A = 0.08 X 10 6 
second. 

1.9 Observed Frequency Variations in Transmittance 

In Fig. 6 is indicated the shapes of the envelope vs frequency varia- 
tions that can be obtained from (3) when the components U and V 
are given by (5) and (6). These fluctuations will vary with time but 
will have the characteristic shapes indicated in Fig. 6, which resemble 
shapes obtained in sweep-frequency measurements on a link of the 
length for which the above value of A applies. 

A better indication of the adequacy of the present idealized tropo- 
scatter model is obtained by comparing the autocorrelation coefficient 
of the envelope as given by (24) with the correlation coefficient derived 
from observations. In Fig. 7 is shown the theoretical coefficient for 
A = 0.08 X 10~ 6 second together with coefficients obtained from three 
experimental runs considered representative. 

The bandwidth capability can be defined as the maximum baseband 
signal spectrum that can be received with some coherence between 
spectral components at the maximum and minimum frequencies. This 




-277/A -77/A 



77/A 



277/A 



Fig. 6 — Illustrative rectified envelope vs frequency characteristic r(u) ob- 
tained with expressions (5) and (6) in (3). The amplitudes Cj at the radian fre- 
quencies uj = jrr/A from the carrier are c, = (a,-* + by*)*. The amplitude of the 
envelope at any intermediate frequency u depends on the amplitudes and phases 
of all cj between j — — oo and j = °o . In sweep-frequency measurements with a 
radian frequency sweep from - tt/A to tt/A from the carrier, the envelope varia- 
tions might be like that in any of the intervals a-b, b-c, c-d, etc. 
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Fig. 7 — Theoretical vs observed envelope autocorrelation functions. Above: 
autocorrelation coefficient obtained from (24) with A = 0.08 X 10 -6 second. Be- 
low: autocorrelation coefficients given in Fig. 70 of Ref. 2 and derived from meas- 
urements of envelope variations with narrow-beam antennas on four days: 1. 
Sept. 13, 1957; 2. Sept. 30, 1957 (considered very unusual); 3. Oct. 15, 1957, and 
4. Nov 8, 1957. The value of A derived from (25) for the experimental link is A = 
0.08 X 10~ 6 second. 



bandwidth is equal to the separation between Cj and c,-+i in Fig. 6, which 
corresponds to the separation between null points in (23), for which 
k(v) = and k t (v) = ir/4. It is given by 1/2A cps and for A = 
0.08 X 10 -6 second is 6.3 mc/second. 

With a smaller spectral bandwidth, distortion will be reduced and 
transmission performance improved. A more realistic appraisal might 
be half the above maximum bandwidth, or 3.15 mc/second, for which 
\ r (f) = 0.9. In Kef. 2 the criterion k 2 ( v) = 0.6 corresponding to K r ( v) = 
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0.904 has been selected, and twice this spectrum bandwidth as required 
in double sideband transmission is quoted in Table VII of the reference. 
The mathematical model represented by (3) to (6) is based on certain 
idealizations outlined in Section 1.7 and in the Appendix. It appears 
from the above that certain theoretical transmittance variations based 
on this model conform sufficiently well with observed variations for the 
model to be acceptable.* In order to determine expected performance 
with digital transmission, it is necessary to consider certain other sta- 
tistical properties of tropospheric channels based on the above model, 
as discussed in sections that follow. 

II. TRANSMITTANCE VARIATIONS WITH TIME 

2.1 General 

As discussed in Section 1.2, the transmission vs frequency charac- 
teristic of a tropospheric scatter channel is a highly variable quantity, 
as indicated in Fig. 1. One way of avoiding transmission impairments 
owing to variations in transmittance with frequency is to transmit by 
narrow-band modulation of a number of different carriers. The amplitude 
vs frequency characteristic can then be regarded as virtually constant 
over each narrow band, and the phase characteristic as linear, as indi- 
cated in Fig. 1. With this method, it is permissible to assume flat fading 
within each narrow band, but the various narrow channels will not fade 
independently. In addition to such flat fading there will be variations in 
the phase and frequency of each received carrier with time. Owing to 
the narrow bandwidth of each channel, the duration T of a signal or 
sampling interval may be relatively long, and it becomes necessary to 
consider the above amplitude, phase and frequency variations over this 
interval T. The probability distribution of these variations are basic 
to later considerations of various digital transmission methods and are 
discussed here. They can be obtained from expressions given by Rice 
for narrow-band random noise. 12 

2.2 Amplitude and Phase Distributions 

Let the frequency w and thus u = co — w be fixed, and consider only 
time variations in r and <p. The probability density of (p is simply p(<p) = 
1/2t, since each phase is equally probable. Since the components U and 
V are the sum of a very large number of independent random variables, 
in accordance with (5) and (6), each component U and V will have a 

* This conclusion appears to be supported by the results of recent measure- 
ments of k(») for a 100-mile path. 24 
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normal law or Gaussian probability density. The probability density 
of the envelope in this case follows the Rayleigh law, and the probability 
that the envelope r exceeds a specified value 1\ is given by 

P(r ^ n) = exp(-r,7f 2 ) (27) 

where r is the rins amplitude of the envelope or the transmittance taken 
over an appropriately long time. 

The average received envelope power is in this case f 2 = S = 2S, 
where S is the average carrier power, i.e., the average power within 
the envelope. The probability that the received envelope power at any 
instant exceeds a specified value <Sj = 2»Si is 

P(S > Si) = exp(-&/<S) = exp(-,S,/.S*). (28) 

The median value S m of »S is obtained from P(S ^ S m ) = \, which 
gives S m — S In 2. Hence, in terms of the median value 

P(S ^ Si) = exp[-WS m ) In 2]. (29) 

The distribution represented by (28) or (29) is shown in Fig. 8. 

The above distribution of rapid fades is to be distinguished from the 
distribution of slow variations in the envelope, or in attenuation, dis- 
cussed in Section 1 .4. 



2.3 Distribution of Envelope Slopes (/•') 

One measure of the rapidity of the above amplitude variations is the 
fading bandwidth discussed in Section l.G. From this fading bandwidth 
can be derived the probability distribution of the slope r' = clr(t)/dl 
in the envelope. 

The rapidity of changes in the envelope and phase depends on the 
time rate of change in the heterogeneities in the common volume — that 
is to say, the variations with respect to time of the coefficients aj(t) 
and bj(t) in (5) and (6). These changes are characterized by the auto- 
correlation function of U(t) and V(t), or by the corresponding power 
spectrum. When the power spectra of U and V are the same, and are 
specified, the probability distribution of /•' = clr(t)/dt and <p' = d<p(t)/dl 
can be determined. These distributions are the same as for random noise 
of specified power spectrum. The probability that | r' \ exceeds a speci- 
fied value | r\ \ follows the normal law " 

P( \r'\ ^ |n' | ) = erfc(fe/2*) (30) 



1(52 
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Fig. 8 — Rayleigh probability distribution of rapid fluctuations in envelope of 
a received carrier owing to multipath propagation. 



in which 



where 
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k = n'/f 

f' = rms amplitude of r' 

= [Uh ~ &i7M* (31) 

bn = f W(y)y n dy. C32) 



The above result (30) follows from equation (4.6) in lief. 12 for Q = 0, 
by integration with respect to 72 = r between and <», and in turn 
with respect to W = r' between r/ and <* . 
Expression (30) can alternatively be written 

P[\r'\ ^ At'] = erfc (fc/2*). (33) 

In the particular case of flat power spectrum W(y) = W of band- 
width y, (32) gives 

/ >0 = Wi; lh = W-f/2- b 2 = Wtf/3 

and (31) becomes 

f' = ry/6* « 0.405ry. (34) 

The fading bandwidth in the above case is y radians/second. 
With a Gaussian spectrum (17) expression (32) gives 

b = tf(0); In = a(2/V)**(0); b 2 = <r 2 *(0) 



1 l^ 5 



and (31 ) becomes 



« 0.42fcr « 0.34f7 (35) 

where y is the equivalent bandwidth given by (19). 

2.4 Distribution of Phase Derivative (<p') 

In considering a small phase change A<p, and over a small interval 
At, it is legitimate to use the probability distribution of the phase 
derivative <p' = A<p/Ar, which is given by [Section 5 of Ref. 12] 

P(|„'| fc|*'|)= 1- ^=y 2 (36) 
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in which 



k = (60/62) V = (6o/6 2 r(A^/Ar) (37) 



where 6 and 6» are given by (32). 

Expression (36) can alternatively be written 

k 

P(| v 1 ^ Hh/b y) = i - ^yt~v- 

«S for fc»l. 



(38) 



For a flat power spectrum W(y) = W of bandwidth y 

(6 2 /6 ) 5 = t73 s ^ 0.58-y. (39) 

For a Gaussian spectrum (17) 

(62/60)' = cr^O.87 (40) 

where 7 is the equivalent bandwidth given by (19). 

2.5 Distribution of Frequency Derivative (<p") 

The probability of exceeding a small variation Ao> in frequency over 
a brief interval At can be determined from the probability distribution 
of*>" = Aw/Ar. 

The probability that <p" exceeds a specified value <p\ is given by 

P(\<p"\ ^ \W) = P(\<P"\ ^ MM 

_2k r dx 

' " t io [g(x) + k*\g(x) (41) 



--[ 



1 tan~ l (k/g\x)) 
(1 + x 2 )i 



dx 



where 



k = b m "/b, (42) 

g(x) = (a - 1 + 4.x 2 ) (1 + x 2 ) (43) 

a = 6064/62 2 . (44) 

Expression (41) is obtained from relation (6.10) of Ref. 12 for 
p(r,<p,<p',<p") for Q = 0, by integration with respect to r, <p and <p', be- 
tween and oc , and 2?r and - <* and + co , respectively, and in turn 
by integration with respect to <p" between <pi" and =° . Considerable 
simplification is required to obtain (41). 
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For very large values of k the following approximation applies 

\ + in (I 



P(\<p"\ ^ kb 2 /h) «-^ 



(i-)] 



where In = log e . 

For a flat spectrum W(y) = W of bandwidth y 

a = 9/5 and b 2 /b = y/Z. 

For a Gaussian power spectrum (18) 

a = 3 and b^/bo = a . 
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(45) 



(46) 



(47) 



The quantity (b 2 /6o)* is the rms frequency of the power spectrum 
and b 2 /b is the "variance." 

The probability distribution (41 ) as obtained by numerical integra- 
tion is shown in Tables II and III for flat and Gaussian power spectra. 
For large values of k, approximation (45) is shown in parentheses. 
These probability distributions are shown in Fig. 9. 

III. TRAN'SMITTANCE VARIATIONS WITH FREQUENCY 



3.1 General 

In the previous section a sufficiently narrow signal band spectrum 
was assumed such that amplitude and phase distortion over the narrow 
band could be neglected. In this case it was necessary to consider time 
fluctuations in the transmittance over a pulse duration T that would be 
relatively long owing to the narrow spectrum bandwidth. 

The other extreme of wideband transmission will now be considered, 
in which the duration of a pulse would be short enough for fluctuations 
in transmittance over a pulse interval to be disregarded. In this case 
it becomes necessary to consider variations in the transmittance with 
frequency over the much greater signal spectrum band. The variations 
in the amplitude and phase characteristics with frequency will fluctuate 
with time, so that it becomes necessary to determine the resultant 

Table II — Probability Distribution P( \ <p" | > ky-/?>) 
for Flat Power Spectrum 



k = o 
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.269 


.238 


.158 


.100 


.051 


.031 (.03) 
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Table III — Probability Distribution P( \ <p" \ > ha 3 ) 
for Gaussian Power Spectrum 



k = 


1 
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3 
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20 


50 


100 


1 


.595 


.447 


.369 


.317 


.280 


.182 


.113 


.057 


.033 (.03) 



transmission impairments on the basis of certain probability distribu- 
tions. 

In a first approximation the departure from a constant amplitude vs 
frequency characteristic will be a characteristic with a linear slope, as 
indicated in Fig. 10, that will vary with time. Similarly the departure 
from a constant transmission delay over the channel band can be approxi- 
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Fig. 9 — Probability that <p exceeds "variance" of fading power spectrum by 
factor k for flat power spectrum with bandwidth y and "variance" -y 2 /3 and for 
Gaussian power spectrum with "variance" a 2 . 
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Fig. 10 — First approximations to random departures from constant amplitude 
and delay characteristics are represented by linear variations with frequency. 

mated by a linear variation in transmission delay. The probability 
distributions of the slopes of these linear variations in the amplitude 
and delay characteristics are the same as for corresponding variations 
with time, with appropriate modification of the basic parameters, as 

discussed in the following. 

3.2 Amplitude and Phase Distributions 

Let the time t be fixed, and consider only variations in r and <p with 
the frequency w of a number of transmitted sine waves. 

Each sine wave could be regarded as a spectral component of a carrier 
pulse of very short duration with an essentially flat and continuous 
spectrum about the carrier frequency. In this case u rather than t is 
changed in expressions (5) and (6) for the two components U(u,t) and 
V(u,t). There will in this case be a particular variation with u for each 
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time t. When observations are made for a sufficiently large number 
of specified times, the resultant probability distribution of the amplitude 
and phase will be the same as discussed in Section 2.2 for variation in 
time for a given frequency u. 

3.3 Slope in Amplitude Characteristic (r) 

At a particular time, the envelopes r(u,t) of the received sine waves 
will vary with frequency u. The slope of the envelope will be designated 
dr(u,t)/du = f. It will have a probability distribution as given by (30) 
for the time rate of change in r(u,t). This probability distribution is 

P( | f | > | h | ) = P( | r | ^ kf) = erfc (k/2 l ) (48) 

where erfc is the error function complement and 

k = h/f 

f = mis value of /• 

= [K&s - fciVM 5 (49) 

except that now 

6„ = ! W{h)& n db (50) 

Jo 

where W(8) is the power spectrum given by (21). When W(8) is given 
by (22), (50) gives 6„ = *(0)/A; 6i = *(0)A 2 /2; 6 3 = *(0)A 3 /3 and 
(49) yields 

f = rA/6* (51) 

where /" = ^(0) 5 is the rms amplitude of the envelope. 

3.4 Envelope Delay Distribution 

The envelope delay at a particular time t and frequency u is given 
by <p = d(p(u,t)/du. The probability distribution of this delay (p is given 
by (36) or (38). Thus 

P(\<p\> \<Pi\) = P\\<p\^ Kbo/b )>] 

k (52) 

" J Vl + /c 2 
where as before 

k = (h/btfyi (53) 

where b and b 2 are given by (50). 
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For a flat power spectrum (22) 

(62/60)* = A/3*«0.58A. (54) 

3.5 Distribution of Linear Delay Distortion 

The slope ip = d<p/dit at a particular time represents linear delay 
distortion. The probability that ip exceeds a specified value ip, is given 
by (41), or 

P(\v\ > l£i|) = Fi\v\ ^ M 2 /M 

_ 2A; r d» 



tt Jo fo(*) + fc*)^*) (55) 

2 r tfuT 1 (k/g*(x)) , 



2 r tan" 1 

TT Jo (1 



dx. 



(1 + *»)» 
For very large values of k (45) applies, or 

1 + In Q + 1 Yl (5G) 



P(|£| ^ fc& 2 /6 ) ^4 



where now 



tt/; 



fc = 60^1/62 (57) 

g ( x ) = (a - 1 + 4.r 2 )(l + .r 2 ) (58) 

a = bjbi/bi (59) 

and 6„ is given by (50). 

For a flat power spectrum ( 22 ) 

b,/b = A 2 /3. (60) 

The probability distribution (55) as a function of h is given previously 
in Table II for a flat power spectrum and is shown in Fig. 9. 

IV. ERRORS FROM TRAXSMITTANCE VARIATIONS WITH FREQUENCY 

4.1 ( lateral 

As discussed later, the error probability in digital transmission over 
noisy channels with selective Rayleigh fading can be approximated by 
combining the probability of errors from three basic sources. One of 
these is errors from random noise determined in the presence of flat 
Rayleigh fading. The second source is errors from time variations in the 
transmittance, which is important at low transmission rates. The third 
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source is errors from transmittance variations with frequency, which 
becomes important at high transmission rates and puts an upper bound 
on the transmission rate for a specified error probability. In this section 
an approximate evaluation is made of errors on the latter account. 
As a first approximation, the statistical properties of transmittance 
variations with frequency, ordinarily referred to as selective fading, can 
be represented by the probability distribution (48) of rand (55) of (£). 
The first of these represents a linear slope on the amplitude vs frequency 
characteristics, and the second represents a linear variation in trans- 
mission delay. Errors will occur even in the absence of noise, when r 
or <f> exceeds certain maximum values. These maxima will depend on 
the spectrum of pulses in the absence of distortion, on the pattern of 
transmitted pulses and on the carrier modulation method. After these 
maximum values are determined, it is possible to determine the proba- 
bility of encountering them with the aid of the probability distributions 
of r and <p given in Section III. 

4.2 Carrier Pulse Transmission Characteristics 

It will be assumed that a carrier pulse of rectangular or other suitable 
envelope is applied at the transmitting end of a bandpass channel. The 
received pulse with carrier frequency m can then be written in the 
general form ' 

P (t) = cos (coot - $ )Ro(t) + sin (wo* - *o)Qo(0 (61) 

= cos M - h- Mt)]Po(t), (62) 

where 

P (t) = [Ro(t) + Qo(t)]\ (63) 

*o(0 = tan -1 [Qo(t)/R«(t)], (64) 

Ro(t) = Po(t) 0OB*(f)> ( 65 ) 

Qo(0 = Po(t) Binft(f). ( 66 ) 

In the above relations Ro and Q are the in-phase and quadrature 
components of the received carrier pulse and P (t) the resultant enve- 
lope. The time t is taken with respect to a conveniently chosen origin, 
for example the midpoint of a pulse interval or the instant at which 
Ro(t) or Po(t) reaches a maximum value. 

Let So(u) be the spectrum of received pulses at the output of the 
receiving filter, i.e., at the detector input, and \f/ (u) the phase function 
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of the spectrum, as illustrated in Fig. 11. The functions R (t) and 
Qo(t) are then given by 13 

Ro = R ~ + Ro + , Qo = Qo~ - Q + , 
Ro = - S ( -u) cos [ut + *„ (-u)] du, (67) 

Ro + = - So(u) cos [ut - Vo(u)] du, (68) 

T >>Q 

i r u o 
Qo~ = - S ( -u) sin [ut + ^ ( -u)] du, (69) 

T J (l 
1 f°° 

Qo = - So(u) sin[ut - V (u)]du. (70) 

The upper limit o> can ordinarily be replaced by oo ? since S ( —wo) = 0. 
Let S(u) be the spectrum in the absence of amplitude distortion, 
and A(u) the amplitude characteristic of the channel. The received 
spectrum is then, for a time invariant channel 

S (u) = S(u)A(u). (71) 

4.3 Ideal Pulse Spectra and Pulse Shapes 

In carrier pulse transmission over an ideal channel the sideband 
spectrum of carrier pulses at the detector input will be symmetrical 
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Fig. 11 — Amplitude and phase functions of pulse spectrum at channel output, 
i.e., detector input. 
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about the carrier frequency. As discussed elsewhere, 1 it is possible to 
realize optimum performance in binary transmission by AM, PM and 
FM with an infinite variety of pulse spectra at the detector input, with 
the general properties illustrated in Fig. 12. With all of these spectra, 
pulses can be transmitted without intersymbol interference at intervals 

T = x/fi = 1/2B (72) 

where B is the mean bandwidth in cps to each side of the carrier fre- 
quency, as indicated in Fig. 12. 

A desirable pulse spectrum in various respects is a raised cosine 
spectrum as illustrated in Fig. 13, given by 



««>-«-)- §«*H- 



(73) 




Fig. 12 — General properties of ideal spectra of carrier pulses at channel out- 
put (detector input) that permit pulse transmission without intersymbol inter- 
ference at intervals T = ir/il = \/2B. 
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Fig. 13 — (a) Raised cosine bandpass pulse spectrum and (1>) carrier pulse 
transmission characteristic, i.e., envelope of a carrier pulse. 



The corresponding carrier pulse at the detector input as shown in 
Fig. 13 is given by 

Po(0 = A>(0 cos (out - <Pq) 

where 

sin fit cos fit 



PM = »o(0 = 



fit 1 - (flt/ir) 2 



(74) 
(75) 



4.4 Linear Variation in Amplitude Characteristic 
Let \po(u) = and 

A(ii) = 1 + c.u 
where c is a constant. In this case ( 71 ) becomes 
S (u) = S(w)(l + cm). 



(76) 
(77) 
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When the received spectrum in the absence of distortion has even 
symmetry about the carrier frequency w , such that S(— u) = S(u), 
(77) in (67) to (70) gives 

R (t) = - f S(u) cos ut du (78) 

7T •'0 

Q (t) = -- f uS(u) sm utdu (79) 

■K •'0 

= c^Ro(t) = cRo(t). (80) 

at 

In the case of a raised cosine spectrum, R (i) is given by (75) and 
(80) yields 

A , A rt ^ cos2fi£ on sin 2SU , * 

ft (0 = c212 2^[i - (2«/>)i " c2 ° (2Q0'[i - BBggff (81) 

= for i = 0. (82) 

At the first sampling points before and after t = 0, t = ±T = 
±(r/G) and(81) yields 

Qo(±7 T ) = ±c«/3t. (83) 

At the next sampling points t = ±2T = ±27r/& 

Qo(±2T) = ±cB/30ir. (84) 

From (83) and (84) it appears that only the first sampling points 
I = ±T need to be considered in determining the effect of linear am- 
plitude distortion. 

4.5 Probability of Errors from Linear Amplitude Distortion 
The rms amplitude of the component Q a (±T) is given by 

Q (±T) = 80/3* = cB/3 (85) 

where B = 2fl/27r and c is the rms amplitude of r as given by (51) or 

c = f = fA/6*. (86) 

Thus (85) becomes 

Qo(±T) = f(£A/3-6 J ). (87) 
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The rms amplitude of 7?o(0) is f. Hence 

OoW BA 
V 4(0) 3-6*- (88) 

This is the ratio of rms intersymbol interference at the first sampling 
points to the rms value of the peak pulse amplitude. 

The probability of exceeding the above ratio by a factor k is, in ac- 
cordance with (48) 

P(u ^ kij) = erfc (k/2 l ). (89) 

The probability of error will depend on the carrier modulation method. 
In general, however, the approximate allowable peak value of rj in the 
absence of noise is 

3 « 5- (90) 

The probability of exceeding this value, corresponding to 
k = 3 • 672dA n 

P € = erfc (3-3*/26a) « erfc (2.6/fiA). (91) 

This probability is much smaller than that resulting from a linear 
variation in delay over the transmission band. For example, if B = 10 6 
cps and A = 10~ 7 sec, 1/BA = 10" 1 and P e = erfc (20), which is negligi- 
ble. 

4.6 Linear Variation in Envelope Delay 

It will be assumed that the phase distortion component is given by 

¥o(u) = cu, (92) 

which corresponds to a linear delay distortion given by 

tfo'(u) = 2cm. (93) 

In this case expressions (07) to (70) give for a raised cosine spectrum 

Ro( — t) = Bo(0 = - / cos x cos ax cos bx 2 dx (94) 

7T ''O 

4 r /2 2 
Qo(—l) = Qo(0 = - I cos x cos «.r. sin bx 2 dx, (95) 

where 

a = 4(7/7'), ^ - (4/*) (tf/7 1 ); 7' = (1/6) 

in which the delay r2 is defined as in Fig. 14. 
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-PULSE SPECTRUM 




LINEAR DELAY 
DISTORTION 

Fig. 14 — Raised cosine pulse spectrum with linear delay distortion. 

The above integrals have been evaluated by numerical integration 
and are tabulated elsewhere. 13 The functions Ro(t) and Q (t) are shown 
in Fig. 15, as a function of t/T = IB for various values of d/T = dB. 
The phase has been adjusted to at t = 0, hence the notation #00 and 

Von • 

4.7 Maximum Tolerable Linear Delay Distortion 

Intersymbol interference at sampling points owing to linear delay 
distortion is significantly greater than that resulting from a linear slope 
in the amplitude characteristic. Moreover, pulse patterns that cause 
maximum intersymbol interference with linear delay distortion will not 
give rise to intersymbol interference from a linear slope in the am- 
plitude characteristic, and conversely. For this reason it suffices to 
consider the more important component, i.e., linear delay distortion. 

The reduction in tolerable noise power owing to linear delay distor- 
tion has been determined elsewhere 13 for binary AM with envelope 
detection, binary PM with synchronous detection, and binary FM with 
frequency discriminator detection. For these methods the reduction 
in noise margin is shown in Fig. 16 as a function of the parameter 
X = d/T = dB. In the same figure is shown the reduction in noise 
margin for two-phase and four-phase modulation, with differential 
phase detection as determined by methods similar to those for the other 
modulation methods in the above reference. These methods essentially 
consist in determining the maximum intersymbol interference that can 
be encountered, considering the pulse shapes shown in Fig. 15 and all 
possible pulse patterns over the number of pulse intervals that contribute 
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4.0 



Fig. 15 — Carrier pulse transmission characteristics for raised cosine pulse 
spectrum and linear delay distortion. For negative values of t/T = IB the char- 
acteristics are the same as shown for positive values. 

significantly to intersymbol interference. Exact analytic determination 
of the maximum impairments does not appear feasible, and it becomes 
necessary to resort to trials for selection of the worst condition. It should 
be noted that with binary PM with differential phase detection the 
optimum threshold level differs from zero owing to a bias component 
in the demodulator output. 13 The curve in Fig. 16 and the analysis that 
follows assume automatic adjustment to the optimum threshold level, 
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1.8 2.0 2.2 2.4 2.6 



Fig. 16 — Maximum reduction in noise margin owing to linear delay distor- 
tion: 1, binary AM with envelope detection; 2, binary FM with frequency dis- 
criminator detection; 3, binary PM with differential phase detection; 4, binary 
PM with synchronous detection; 5, four-phase modulation with synchronous 
detection; 6, four-phase modulation with differential phase detection. 

and a significantly greater error probability would be encountered with 
zero threshold level. 

It will be noted that the noise margin is reduced to zero for certain 
values X of X. These values apply for certain combinations of baseband 
pulses in about four pulse positions. The probability of this and other 
pulse patterns must be considered in evaluating error probability as 
discussed below. 

4.8 Probability of Errors from Linear Delay Distortion 

As X is increased slightly above the value Xo mentioned above, inter- 
symbol interference increases rapidly. Thus errors will occur for a value 
X, of X only slightly greater than Xo , for certain combinations of two 



I)I< ilTAL TKOPO.SCATTER 



17'.) 



pulses, occurring at times —T and +7' relative to the sampling instant 
I = 0. There are four possible combinations of these two pulses. For 
one of these (say 1,1), an error will occur if X ^ X,. . For another (say 
— 1, —1), an error will occur if X ^ — X, . For the other combinations 
( — 1, 1) and (1, —1), intersymhol interference will cancel so that the 
probability of error is zero. The probability of error is thus 



1\ = i(i + \)P( |X| ^ |X e | ) 

= \P( i x : ^ i x, i > 



(96) 



where P( | X | ^ | X,- | ) is the probability that the absolute value of 
X is greater than \ r . 

For a given value X, = d,B the corresponding slope <p of the linear 
delay distortion is 



!p e = d r /2irB 

= X„/2tB 2 . 
The following relation applies 

P( I X I ^ | X, I ) = P( I ip I ^ 



). 



(97) 



(98) 



The probability distribution represented by the right-hand side of 
(1)8) is given by (55) with v?i = ip t: . For small probabilities (56) applies, 
so that in view of (96) and (08) the error probability is 



i r 

= 2*k e _ 



<Pr 



) 



K 



1 --ln(- e + 1 



where 



= 3X, 2irXB' 



Witli (100) in (00 



/', = 



fh'-O 



3X r 



1tA-B-/_ 



(00) 



( 100) 



(101) 



From Fig. 10 it will lie noted that for binary AM and FM, and for 
binary I'M with differential phase detection, X ~ 1.8. For these cases 
it appears a legitimate approximation to take \ r = 2. On this premise 
the error probabilities given in Table IV are obtained for various values 
of the parameter <\B. 
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Table IV — Probability of Errors in a Digit Owing to Linear 

Delay Distortion in Absence of Noise for Binary AM, FM 

and PM (with Differential Phase Detection) 



ab = lir* 


10-3 


10-2 


io-' 


3.1 X 10- 8 


2.4 X 10-° 


1.6 X 10~' 


8 X IO" 3 



The above error probabilities are shown in Fig. 17 as a function of 
AB. If, for example, A = 10~ 7 second and B = 10 5 cps, then AB = 10~ 2 
and P r = 1.6 X 10 -4 . Pulses could in this case be transmitted at a 
rate of 100,000 per second with a minimum error probability P e = 1.6 X 
10~\ In the presence of noise the error probability will be greater, as 
discussed in a later section. 
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Fig. 17 — Error probability in binary AM, FM and PM owing to linear delay 
distortion for maximum departure A (seconds) from mean transmission delay. 
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The error probability with four-phase modulation and differential 
phase detection can be determined in a similar way. In this case A ~ 0.9 
and A,^ 1 in (101). 

V. ERRORS FROM TRANSMITTANCE VARIATIONS WITH TIME 

5.1 General 

As mentioned in Section 4.1, transmittance variations with time is a 
second basic source of error in digital transmission. In transmission at 
low rates the bandwidth B of the pulse spectra will be narrow, so that 
fading can be regarded as constant over the spectrum band. Errors 
from selective fading, as considered in Section IV, can then be disre- 
garded. On the other hand, the duration of a signal interval T may then 
be sufficiently long so that consideration must be given to random 
fluctuations in the amplitude, phase and frequency of the carrier between 
one signal interval and the next. Errors may occur owing to such fluctua- 
tions even in the absence of noise. The probability of errors in this 
account is evaluated here. 

5.2 Amplitude Variations 

The amplitude of a received wave will fluctuate with a Rayleigh 
distribution (10). Because of the great range of fluctuation, it is essential 
to provide automatic gain control at the receiver to prevent overloading 
and resultant adverse effects. Such gain control is activated by circuitry 
that integrates the received wave over a number of signal intervals T. 
With FM and PM only a few pulse intervals are required, for the reason 
that the received carrier wave is essentially independent of the pulse 
patterns. It is thus possible to provide effective gain control against rapid 
variations in the received carrier wave that occurs over a few signal 
intervals. Moreover, with FM and PM the distinction between marks 
and spaces is made by positive and negative deviations from zero thresh- 
old level in the detection process. This permits the use of limiters at 
the input to the detectors, to prevent the adverse effect of rapid fluctua- 
tions in the amplitude of the received carrier wave owing to fading. These 
advantages in applications to fading channels are not shared by AM, 
for reasons outlined below. 

In binary AM or on-off carrier transmission, the received wave may 
be absent over a large number of consecutive signal intervals T. Hence 
automatic gain control must be activated by circuitry that integrates 
the received pulse train over a very large number of signal intervals T; 
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otherwise gain would be increased during long spaces, regardless of the 
fading condition. For this reason automatic gain control is inherently 
slow, in relation to the duration of a signal interval. It may thus be 
ineffective as applied to transmission at slow rates. With transmission 
at high rates, however, such that variations in the received wave owing 
to fading are inappreciable even over a large number of signal intervals, 
it may be possible to implement effective gain control. 

At low transmission rates, such that fading is virtually constant over 
the band of the pulse spectrum, intersymbol interference can be made 
inappreciable. In this case it is possible to employ limiting prior to 
detection, and this method may then be more effective than automatic 
gain control, or could be used in conjunction with it. The limiter would 
slice the received wave at an appropriately selected level L. In the 
choice of the optimum slicing level it is necessary to consider the proba- 
bility of errors during a mark owing to fading such that the received 
wave is less than L. In accordance with (10") this probability is 

P(r^L) = 1 -exp(-L 2 /f 2 ) 
« L 2 /f\ 

A second consideration in the choice of L is the probability of errors 
owing to noise during a space, which is increased as L is reduced. The 
optimum threshold level considering both effects is determined in Sec- 
tion 6.9. 

Owing to even small intersymbol interference, the use of a limiter as 
postulated above may be precluded in actual systems. For example, 
let L be 10 per cent of the rms signal amplitude f, and let intersymbol 
interference be 5 per cent of L when the received signal is just equal to 
L. When the received signal is increased by a factor 20, intersymbol 
interference would be increased correspondingly and would be equal to 
L. Hence errors would occur even in the absence of noise. This is the 
inherent reason why limiting is generally ineffective as applied to binary 
AM. However, even if intersymbol interference could he disregarded, 
the error probability in the presence of noise will be greater than with 
binary PM or FM, as shown in Section 6.9. 

5.3 Carrier Frequency Variations 

In transmission over troposcatter links, random fluctuations will 
occur in the carrier frequency, which may be important from the stand- 
point of receiver implementation with any modulation method. Such 
fluctuations can be limited at the input to the IF filter with the aid of 



DIGITAL TROPOSCATTER 183 

signal-tracking oscillators for demodulation of the received radio fre- 
quency wave. The frequency of such oscillators may be controlled by 
feedback from the mixer output or from the detector output. The follow- 
ing expressions apply for the probability distribution of carrier fre- 
quency fluctuations without such frequency control at the receiver. 

The probability distribution of frequency variations is given by 
(38). For a Gaussian fading power spectrum, the probability that the 
frequency variation <p' = Aco exceeds k<r is thus 

P( | Aw | £ hr) « (1/2/c 2 ). (103) 

The equivalent fading bandwidth is in accordance with (19) y ~ 1.25c - . 
The probability that Aco exceeds ky is thus 

P( | Aco | ^ ky) » (I /3k 2 ). (104) 

Since a and y are nearly proportional to the carrier frequency, it 
follows that the frequency fluctuations encountered with a specified 
probability will be nearly proportional to the carrier frequency. By 
way of example let y ~ 2 radians/second or about 0.3 cps. The proba- 
bility that the frequency fluctuation exceeds 30 cps is in this case ob- 
tained from (104) with k = 100 and is 3 X 10 -5 . It appears that for 
bandwidths of the pulse spectra in excess of about 5000 cps, frequency 
fluctuations will not be important. However, for narrow band spectra 
the random frequency excursions may become excessive and give rise 
to errors, particularly with frequency modulation, as discussed below. 

5.4 Frequency Variations over a Signal Interval 

It will be assumed that the carrier frequency excursion is limited with 
the aid of a signal-tracking oscillator, or that a demodulation process is 
used in binary FM in which the change from mark to space is based on 
comparison of the frequencies in adjacent signal intervals of duration 
T. If the separation between mark and space frequencies is 2O i , an 
error will occur if the frequency is changed by +fioi for a space and 
by — ftoi for a mark. 

From (41 ) it is possible to determine the probability of errors owing 
to frequency changes ±S2oi over a signal interval of duration T. The 
maximum permissible value of <p" is determined from 

¥W*T = ±Om (105) 

where the positive sign applies for a space and the negative sign for a 
mark. 
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With an ideal pulse spectrum the pulse interval is given by T = r/O, 
so that (105) can be written 



<f" max = ifioi^/'T- 



(106) 



5.5 Error Probability in Binary FM 
The error probability is in this case 

P t = \P{ \ V »\> |^nax" |) (107) 

where the factor \ occurs when the probability functions is defined in 
terms of the absolute values as in (41). 

The parameter k defined by (42) in this case becomes 

ft max — ^max / & 

, 2 ( 108 ) 

= S2oi^/7TO' . 

With frequency discriminator detection, ft i = ft- For a raised cosine 
spectrum, B = 2B = Q/v and 

= *B*/ct\ 



fc, 



(109) 



Employing (45), the probability (107) of an error becomes 

p -=Ci)T i+in ( i+ S)]- (uo) 

In the above relation, a is in radians/second while B is in cps. The 
equivalent fading bandwidth is, in accordance with (19), y W 1.25a. 
The ratio of the maximum bandwidth B in cps to y in cps is thus 



B 



M = - 



2wB 5B 



7/2tt 1.25a 



(111) 



The probability of error (110) is given in Table V for various ratios ft. 
These error probabilities are shown in Fig. 18. 

Table V — Error Probabilities with Binary FM from 
Flat Rayleigh Fading in Absence of Noise 



H = 10 


100 


1000 


10000 


B/a = 2 


20 


200 


2000 
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1.8 X 10-" 
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Fig. 18 — Error probability in binary FM in absence of noise, owing to fre- 
quency variations over a pulse interval T resulting from flat Rayleigh fading. 



5.G Phase Variations over a Signal Interval 

The probability density of the carrier phase is 1/2t, such that any 
phase may be encountered unless the carrier phase wander is limited 
by phase tracking oscillators in the demodulation process. In a digital 
phase modulation system where appreciable phase wander may be 
expected, the preferable demodulation method is differential phase 
detection. With this method the phase error will be limited to that 
encountered over a signal interval T. 

From (36) it is possible to determine the probability of an error for 
a given maximum tolerable phase change 6 over an interval T. For 
k » 1 the following relation applies 
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P(|/| ^ |« |) = ^5 (H2) 



6 2 T 2 



(113) 



26 2 ' 
With a Gaussian fading power spectrum (40) applies and 

P[\v'\ ^ («')] = (.^T 2 /26 2 ). (114) 

5.7 .Error Probabilities in PM 

With two-phase modulation = ±(ir/2), while with four-phase 
modulation 6 = ±(ir/4). Hence the probability of error with these 
methods as obtained from (114) is, for two-phase modulation 

P e tt(2/T 2 )a 2 T 2 tt0.2a 2 T 2 (115) 

and for four-phase modulation 

P e « (8/ttVT 2 « 0.82a 2 T 2 . (116) 

These expressions apply provided the signal duration is sufficiently 
short so that the change in phase is small and can be considered linear 
over the interval. More accurate expressions that do not involve this 
assumption have been derived by Voelcker 9 for the error probability. 
Thus, with two-phase modulation the error probability is actually 

P e = |[1 - k(T)\ (117) 

and with four-phase modulation 

P. - 1 - | «(T)[2 - «T)T* tan"' [2 _'%, )]t ("« 

where k(T) = k(t) for t = T, i.e., the autocorrelation function for 
each quadrature component as defined by (15). 

For a Gaussian fading spectrum, k{T) as obtained from (17) is 

k{T) =exp(-(r 8 T 8 /2). (119) 

For <tT « 1 : 

k(T) Rll - <7 2 r 2 /2. (120) 

With the latter approximation in (117) and (118), the error proba- 
bility with two-phase modulation becomes 

P„^i<r 2 r 2 = 0.25er 2 T S (121) 
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and with four-phase modulation 



1 . l\ 2n* 



Pt = vi + J ; ^ 0,82cr ^ 122 ^ 

which are to be compared with (115) and (116), respectively. The 
somewhat greater inaccuracy with two-phase than with four-phase 
modulation comes about since the phase change ± (71-/2) cannot be 
considered small as required for (114) to apply. 

In the above relations T is the interval between phase changes, which 
is related to the bandwidth of the baseband pulse spectrum. With 
idealized spectra of the type shown in Fig. 12, the interval is 

T = 1/2/? (two-phase) (123) 

= 1/4B (four-phase) (124) 

where B is the equivalent mean bandwidth. 

In the particular case of pulses with a raised cosine spectrum, the 
maximum bandwidth is 

B = 2J5 (125) 

so that 

T = 1/B (two-phase) 

(126) 
= 1/2B (four-phase). 

In terms of the above bandwidth the error probabilities (115) and 
(116) are thus the same for both two-phase and four-phase modula- 
tion and are given by 

P e ^0.05(a/B) 2 (127) 

«0.2(er/£) 2 . (128) 

The above relations apply for any number of phases. For this reason 
the capacity of a noiseless channel could be increased indefinitely by 
increasing the number of phases. There will, however, be certain limita- 
tions in this respect owing to intersymbol interference, as in stable 
channels. 

The above error probability is shown in Table VI for various values 
of B/a and n = bB/a, where ju is the ratio defined by (111). It will be 
noted that these error probabilities are somewhat smaller than with 
binary FM as given in Table V. 

The above probabilities of an error in a single digit are shown in Fig. 
19, as a function of /i. 
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Table VI — Error Probabilities with Differential PM 
from Flat Rayleigh Fading in Absence of Noise 
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As noted in Section 1.6, there will be a certain median value of 7 and 
thus a certain median value of n and corresponding median error proba- 
bility. During certain intervals, the error probabilities will be signifi- 
cantly smaller or significantly greater than the median error proba- 
bilities. 
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Fig. 19 — Error probability in binary PM with differential phase detection 
in absence of noise, owing to phase variations over pulse interval T resulting from 
flat Rayleigh fading. 
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VI. ERRORS FROM NOISE WITH FLAT RAYLEIGH FADING 

6.1 General 

As mentioned in Section 4.1, a third basic source of errors in tropo- 
scatter transmission is random noise. The probability of errors from 
noise depends on the modulation and detection methods and on their 
implementation. For optimum performance it is in the first place neces- 
sary to have appropriate pulse spectra such that intersymbol inter- 
ference is avoided in transmission over ideal channels. Moreover, the 
error probability depends on the division of spectrum shaping between 
transmitting and receiving filters. The minimum error probabilities 
with various modulation and detection methods as quoted here are 
based on optimum design in the above and various other respects, such 
as accurate sampling of pulse trains. The probability of errors from 
noise in actual systems will be greater owing to various imperfections in 
implementation. 

6.2 Signal-to-N oise Ratios 

In carrier pulse transmission over an ideal channel, the sideband 
spectrum of the carrier pulses at the detector input will be symmetrical 
about the carrier frequency. As discussed elsewhere, 14 it is possible to 
realize optimum performance in binary transmission by AM, PM and 
FM with an infinite variety of pulse spectra at the detector input with 
the general properties discussed in Section 4.3. 

The error probability in digital transmission over noisy channels is 
ordinarily specified in terms of the average signal-to-noise ratio at the 
input to the receiving filter that ordinarily precedes the detector. This 
signal-to-noise ratio is ordinarily taken as 

P = S/N 

8 = average carrier power at detector input 

N = average noise power in a flat band B = 1/27' at 
input-to-receiving filter. 

When S represents the average signal power in a fading channel, the 
designation p = S/N will be used in place of p. 

The above reference band B is the minimum possible bandwidth in 
baseband pulse transmission without intersymbol interference. The 
minimum possible bandwidth in double sideband transmission as used 
in binary AM, PM and FM is 2B. 

The error probability as related to p will depend on the division of 
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spectrum shaping between transmitting filters and the receiving filter 
at the detector input. With optimum division, the error probability is 
the same as for transmission over a flat band B to each side of the carrier 
frequency. 14 Such a flat channel band is ordinarily assumed or implied 
in theoretical analyses, though not feasible in actual systems. 

6.3 Error Probabilities with Flat Rayleigh Failing 

Let r be the signal amplitude and P e °(r) the error probability of errors 
owing to random noise in transmission over a stable channel with signal 
amplitude r. In the presence of fading, let the probability density of 
various signal amplitudes be p(r). The error probability in transmis- 
sion over fading channels is then 

P e = f P e °(r)p(r) dr. (129) 

With Rayleigh fading the probability density p(r) is the derivative 
of (27) with respect to r x . With r in place of n the probability density is 

p(r) = (2r/f a )exp(-r 2 /f 8 ) (130) 

= (rAS)exp (r 2 /2S) (131) 

where S = f 2 /2 is the average signal power. 

6.4 Binary PM with Synchronous Detection 

In binary PM, marks and spaces are transmitted by phase reversals. 
With ideal coherent or synchronous detection the error probability 
in transmission over a stable channel is 

P a °= ierfc(p/2) J . (132) 

The error probability with Rayleigh fading as obtained from (129) 
is, in this case 7,0 

'■-ft-UO'Ki (133) 

where p = S/N = ratio of average received signal power with Rayleigh 
fading to average noise power as previously defined. 

6.5 Binary PM with Differential Phase Dection 

With binary PM and differential phase detection the error proba- 
bility in transmission over a stable channel is ° 

P.° = he-". (134) 
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The error probability with Rayleigh fading is, in this case 9 

P e = l/2(p + 1). (135) 

6.6 Binary FM with Dual Filter Detection 

With this method two receiving filters are used, centered on the space 
and mark frequencies wi and co 2 , as indicated in Fig. 20, with sufficient 
separation to avoid mutual interference between the space and mark 
channels. Complementary binary amplitude modulation is used at the 
two carrier frequencies, and the two baseband filter outputs are com- 
bined with reversal in the polarity of one. 

The error probability in transmission over stable channels with co- 
herent detection is 



Pe = \ erfc (p 3 /2) 
and with noncoherent detection is b 

P.° = §exp(-p/2). 



136) 



(137) 



4B = 4/T 




(a) 




(b) 

Fig. 20 — Comparison of channel bandwidth requirements in binary FM with 
(a) frequency discriminator detection and (b) dual filter detection. 
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Comparison of (136) with (132) shows that the error probability 
P e with Rayleigh fading is obtained by replacing in (133) p with p/2. 
This yields for coherent detection 

i. (138) 

Comparison of (137) with (134) shows that P e is obtained by re- 
placing in (135) p with p/2, in which case, for noncoherent detection 

P e =l/(p + 2). (139) 

6.7 Binary FM with Frequency Discriminator Detection 

With this method a single receiving filter is used, with space and mark 
frequencies as indicated in Fig. 20. Pulse transmission without inter- 
symbol interference over a channel of the same bandwidth as required 
for double-sideband AM is in this case possible for certain ideal ampli- 
tude and phase characteristics of the channels, as shown elsewhere. 

The error probabilities in the absence of fading depends on the charac- 
teristics of the bandpass channel niters and the post-detection low-pass 
filter, and are difficult to determine exactly. Approximate evaluations 1 
indicate that for a given error probability, about 4 db greater signal- 
to-noise ratio would be required than for binary PM with coherent 
detection, when no post-detection low-pass filter is used. Recent exact 
evaluations by Bennett and Salz, 17 indicate 3 to 4 db increase in the 
required signal-to-noise ratio over a variety of filter shapes. With an 
optimum post-detection low-pass filter, a small improvement may be 
realized, such that about 3 db increase over binary PM with coherent 
detection would be expected. On this basis it appears that the error 
probability will be virtually the same as for binary FM with dual filter 
coherent detection, such that the principal advantage over the latter 
method is a two-fold reduction in bandwidth. 

6.8 Binary AM with Ideal Gain Control 

It will be assumed that the receiver can be implemented with ideal 
automatic gain control, such that the output in the presence of a mark 
would have a fixed level I and in the presence of a space would be zero. 
This condition can be approached at sufficiently high transmission rates, 
such that the received wave prior to gain control changes insignifi- 
cantly over a large number of pulse intervals of duration T. Under this 
condition the fading bandwidth is negligible relative to the bandwidth 
of the baseband pulse spectrum. 
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On the above premise and with ideal coherent (or synchronous) 
detection, the optimum threshold level for decision between marks and 
spaces would be 1/2. The tolerable peak noise amplitude before an error 
occurs would be 1/2, as compared with I for binary PM, resulting in 
6 db reduction in noise margin. On the other hand, the average trans- 
mitter power is 3 db less than with binary PM. Hence this method 
would have a 3 db disadvantage compared to binary PM with synchro- 
nous detection. 

Accordingly, (132) would be replaced by 

P e ° = |erfc( P /4)* (140) 

and (133) would be replaced by 



P =- 
le 2 



1 - 



+ 2 



(141) 



The above relations are the same as (136) and (138) for binary FM 
with dual filter coherent detection, and (141) is virtually the same as 
(135) for binary PM with differential phase detection. Hence binary 
AM offers no advantage in signal-to-noise ratio even at sufficiently 
high transmission rates such that ideal gain control could be imple- 
mented. 

6.9 Binary AM with Optimum Fixed Threshold Detection 

At low transmission rates, such that the received wave can change 
appreciably over a few pulse intervals owing to fading, gain control 
cannot be effectively implemented, as discussed in Section 5.2. Without 
effective gain control, there will be a certain optimum threshold for 
distinction between marks and spaces. This optimum level and the cor- 
responding signal-to-noise ratio is determined here on the premise that 
no gain control is used. This threshold level could be implemented by 
either a predetection or a postdetection limiter. Assume a probability 
\ of a mark being present; in the absence of noise, the probability of 
errors in marks is, in view of ( 102) 

Pe(r S L) = Ml - exp {-I//2S)] (142) 

where L is the threshold level. In the presence of noise the error proba- 
bility will be only slightly greater than (142). 

A second consideration in the choice of L is the probability of errors 
during a space. This error probability is obtained from (137) with 
p = L /N and is 

P e (n ^ L) = h exp (-L 2 /2N) (143) 
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where n is the instantaneous noise amplitude and N the average noise 

id error probability is 

P e = Ml - exp (-m/2) + exp (-p/x/2)] (144) 



power. 

The combined error probability is 



where 

M = L 2 /S; p = S/N. (145) 

The optimum L or p. is obtained from the condition dP e /dp = 0. This 
yields the following relation for the optimum value p 

exp (-mo/2) = pexp (-Pim/2) (146) 



or 



= 21n ? = 4 ' 606 lQ gio P (147) 

p — 1 p — 1 



In practicable systems p » 1, in the order of 100 or more, and po « 1. 
With (147) in (144), the following approximation is obtained for the 
minimum error probability 

P e . mi „^i[^ I + exp(-lnp)]. (148) 

The above error probability is significantly greater than with binary 
PM or FM. The error probability (148) is thus greater than for binary 
FM with dual filter coherent detection by a factor of at least In p. 
For p = 1000 (30 db) this factor is about In p « 7. Hence about 10 
logio 7 « 8.5 db greater average signal power would be required than 
with binary FM. This assumes that excessive intersymbol interference 
is avoided, which may not be feasible for reasons mentioned in Section 
5.2. Since it is evident that binary AM is at a considerable disadvantage 
in signal-to-noise ratio as compared to binary PM and FM, it will not 
be considered further herein. 

6.10 Combined Rayleigh and Slow Log-Normal Fading 

In the previous determination of error probabilities, rapid Rayleigh 
fading was assumed, with a fixed mean signal-to-noise ratio p over the 
interval under consideration. It will now be assumed that in this interval 
there is a slow log-normal variation in path loss and thus in signal-to- 
uoise ratio at the receiver, in conjunction with rapid Rayleigh fading. 

Let P e be the error probability with Rayleigh fading as previously 
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related to the mean signal-to-uoise ratio p = s/ii, where 5 is the rms 
signal amplitude and n the rms noise amplitude. If p(s) is the proba- 
bility density of the rms amplitudes with slow fading, the probability 
of error in an interval during which the rms amplitude exceeds Si is 

I\i = f P e (s)p(s) ds. (149) 

For p » 1, the expression for P e (s) is of the general form 

P e (s)^c/p =q£r a . (150) 

s 2 /n 2 

For binary PM with differential phase detection and for binary PM 
with coherent dual filter detection, c = J. 

The probability density p(§) is given by (12), or in the present 
notation 

p(«) = —j= — exp [-(In 5/«o) 2 /2(7 2 ] (151 ) 

V 2ir «» 

where s is the median rms amplitude and a is the standard deviation 
of the fluctuation in s. 

With (150) and (151) in (144) 

Pe.i = c -J= - I -1- i exp [- (In s/.s ) 2 /2<7 2 ] tfs (152) 
V2tt °" J '*i r/n* s 

= | -4= f 4oxp [-(* In p/p ) 2 /2a 2 ] rfp (153) 

* V 2tt j pi P" 

where p = s 2 /?I 2 on pi = Si 2 /fl 2 . 
Solution of (153) yields the relation 

P e ,i - P,-tj(<t,k) (154) 

where 

k = Pi/po (155) 

and 

r,(a, k) = § exp (2a 2 ) erfc \t7jT [4<t 2 + In k]\ . (156) 

For pi = 0, In k = — oo and erfc ( — » ) = 2. Hence for this case 

r, = exp (2<r 2 ). (157) 
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This is the factor by which the error probability taken over a long inter- 
val is greater than without a log-normal variation in signal-to-noise 
ratio and only rapid Rayleigh fading. 

Instead of modifying the error probability as above, an alternative 
method is to use an equivalent mean signal-to-noise ratio p e that is 
smaller than p by the factor exp( — 2a ). Thus 

p e = pexp(-2<r 2 ). (158) 

When p e , p and a are all expressed in db, expression (158) can alterna- 
tively be written 

Pe.db = pm — oW/8.09. (159) 

For example, with a representative value <r,ib = 8 db, the last term 
in (159) is 7.4 db. Thus the charts in the later Figs. 21 and 22 apply 
when p is taken 7.4 db less than the median signal-to-noise ratios with 
log-normal fading. 

VII. COMBINED ERROR PROBABILITY 

7.1 General 

In Sections IV to VI, three basic sources of errors in digital transmis- 
sion over troposcatter links were discussed, and expressions were given 
for the probability of error from each of these sources in the absence 
of the others. In a first approximation, the error probability considering 
all three sources can be evaluated by taking the sum of the three error 
probabilities. Approximate expressions are given here for the resultant 
error probabilities, together with charts that facilitate determination 
of error probability as a function of the binary pulse transmission rate, 
when the basic system parameters are known. These are the average 
signal-to-noise ratio p, the mean fading bandwidth y, and the maximum 
departure A from the mean transmission delay. The error probability 
for a given transmission rate can be reduced by various means that may 
or may not entail an increase in total transmitter power or bandwidth 
or both. For a given total transmitter power and bandwidth, the most 
effective means to this end is diversity transmission over independently 
fading paths, as discussed briefly herein. 

7.2 Combined Error Probability 

As a first approximation, the error probability is given by 

P e « P e w + P e w + P. m (160) 
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where 

P e a) = probability of errors in the absence of noise owing to inter- 
symbol interference caused by frequency selective Rayleigh 
fading (Section IV) 

Pe 2) = probability of errors in the absence of noise owing to random 
variations in carrier phase or frequency (Section V) 

P e } = probability or error owing to random noise with nonselective 
Rayleigh fading (Section VI). 

As will be evident from the preceding discussion, and from charts 
that follow, P e {1) can be disregarded when P e (2) must be considered, and 
conversely, for error probabilities P e <3) in the range of practical interest. 
Hence in actual applications (160) will take one of the following forms 

P e ttP e (1) + P e (3) (161) 

P e ttP e (2) + P e (3) . (162) 

In addition, there are intermediate cases in which P e ~ P e <3> . 

In an exact determination of the error probability (161) it is neces- 
sary to consider the net effect of random intersymbol interference on 
the probability of errors owing to random noise, and similarly an exact 
determination of the error probability (162) the probability distribution 
of random phase deviations is involved. Intersymbol interference at a 
particular sampling instant may reduce or increase the tolerance to 
noise, and the net effect considering all pulse patterns may be such that 

(161) is a legitimate approximation. Similarly, random fluctuations 
in the slope of the phase characteristic may decrease or increase the 
tolerance to noise at a particular sampling instant, and the net effect 
considering all sampling instants may be such that (162) is a valid 
approximation. This is evidenced by the following exact relation derived 
by Voelcker 9 in place of (162) for binary PM with differential phase 
detection 

P e = [p/(p+ 1)]P, (2) + /V 3) . (163) 

Since p would ordinarily exceed 100 (20 db), it follows that in this case 

(162) is a very good approximation to (163). 

The exact error probability (161) depends on the probability distribu- 
tion of phase distortion in conjunction with the probability distribution 
of intersymbol interference, which involves consideration of all pulse 
patterns. The combined probability distribution, and in turn the exact 
error probability, would be very difficult to determine, and hence the 
inaccuracy involved in (161) cannot readily be assessed. However, if 
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the probability distribution of intersymbol interference were the same 
as that of the reduction in tolerance to noise owing to random fluctua- 
tions in the slope of the phase characteristic, the inaccuracy in (161) 
would be no greater than that indicated by (162) versus (163). In 
most engineering applications, substantially greater inaccuracy would 
be permissible in the estimation of error probability, such that (161) 
and hence (160) can be considered permissible approximations in the 
present context. 

The above expression (160) is applied below to binary PM and FM. 

7.3 Binary PM with Differential Phase Detection 

For binary PM with differential phase detection P e a) is given by 
(101) with X„ = 2 or 

P (i) _ L 

■L r. — 



[ 1 + h ( 1+ SOT)]- ° 64) 



6 

This error probability is given in Table IV as a function of AB. 

The error probability P, (2) is obtained from (117), or approximation 
(121) 

p« m i[i _ K ( T )} (165) 

« 0.25(<r!T) 2 fef 0.06(<r/B) 2 (166) 

« 0.039(y/B)\ (107) 

The error probability P e w is given by (135) or 

P e (3) = l/2(p+ 1). (168) 



(1) D (2) „.„] r> (3) 



7.4 Error Probability Charts for Biliary PM 

In Fig. 21 are shown the error probabilities P e a \ P. (2) and P 
a function of the transmission rate, for a raised cosine spectrum. The 
error probability P e w depends on the maximum deviation A from the 
mean transmission delay, and curves are shown for a number of values 
of A. The probability P. (2) depends on the mean fading bandwidth y, 
and curves applying for several values of y are shown. Finally, the error 
probability P e i3> depends on p, and is shown for a number of different 
values of p. 

By way of illustration, the combined error probability obtained from 
(170) is shown by the dashed line in Fig. 20 for the particular case in 
which A = 10~ 7 second, y = 2 cps and p = 10 4 (40 db). 

The error probability as a function of transmission rate shown by 
this dashed line could apply to a variety of tropospheric scatter links, 
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Fig. 21 — Probabilities of errors in binary PM with differential phase detec- 
tion: 1, curves for various departures from mean delay show error probabilities 
in absence of noise owing to pulse distortion from selective fading; 2, curves for 
various mean fading bandwidths y show error probabilities in absence of noise 
owing to random phase variations caused by flat fading; 3, curves for various 
mean signal-to-noise ratios p show error probabilities owing to noise for flat Ray- 
leigh fading; 4, dashed curve shows approximate combined error probability for 
p = 40 db, A = 10 -7 second, and 7 = 2 cps. 

since A depends on the length of the link and on the antenna beam 
angles. Moreover, p depends on the transmitter power, the length of 
the link, and the antenna beam angles. Hence, given values of A and p 
can be realized for a great variety of conditions. 



7.5 Binary FM with Frequency Discriminator Detection 

With frequency discriminator detection, the minimum required band- 
width for a given pulse transmission rate is the same as for binary PM, 
and half as great as that required with dual filter detection. 



p (2) ... 
1 e — 
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The error probability P e (1> is in a first approximation the same as 
(161) for binary PM with differential phase detection. For the error 
probability P a w , approximation (110) applies, or 

This error probability is given in Table V as a function of B/<r. 

The probability of error owing to noise is, in a first approximation, 
the same as given by (139) for dual filter detection with coherent detec- 
tion, or 

P e (3) ttl/2p. (170) 

7.6 Error Probability Charts for Binary FM 

In Fig. 22 are shown the error probability P e (1) , P e (2) and P e (3) for 
binary FM as a function of the transmission rate. The curves apply for 
a raised cosine pulse spectrum, and the same basic parameters <r, y and 
p as shown in Fig. 21 for binary PM. The error probability for the partic- 
ular set of parameters previously assumed in Section 7.4 is shown by 
the dashed curve. 

Comparison of the curves in Figs. 21 and 22 shows that the error 
probabilities are the same with both methods except at very low trans- 
mission rates. This applies only as a first approximation and with ideal 
implementation of both methods. 

7.7 Diversity Transmission Methods 

In diversity transmission, either space, frequency or time diversity 
can be used. The performance would be the same with these methods, 
and is an optimum when there is no correlation between the diversity 
paths. This entails adequate separation of receiving antennas in space 
diversity, adequate frequency separation in frequency diversity, or 
adequate time intervals between repetition of signals in time diversity. 

With any one of the above three methods, different combining or 
decision procedures can be used at the receiver, as discussed in considera- 
ble detail by Brennan. 17 The optimum method from the standpoint of 
minimum required signal power for a specified error probability is known 
as "maximal ratio combining," in which the gain of the receiver in each 
path is made proportional to the input signal-to-noise ratio. This method 
is difficult to implement, and a simpler but somewhat less efficient 
method is "equal gain combining," in which the various receivers have 
equal gain and the demodulator baseband output are combined linearly. 
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Fig. 22 — Probabilities of errors in binary FM with frequency discriminator 
detection: 1, curves for various departures A from mean delay show error proba- 
bilities in absence of noise owing to pulse distortion from selective fading; 2, 
curves for various mean fading bandwidths y show error probabilities in absence 
of noise owing to random frequency variations caused by flat fading; 3, curves 
for various mean signal-to-noise ratios p at detector input show error probabili- 
ties owing to noise for flat Rayleigh fading; 4, dashed curve shows approximate 
combined error probability for p = 40 db, A = 10~ 7 second and 7 = 2 cps. 



This entails a demodulator in each diversity channel and common gain 
control of the various channels. The need for a demodulator in each 
diversity channel and common gain control is avoided with "selection 
diversity," in which the receiver having the largest signal is selected. 
Though this method is somewhat less efficient than equal gain combin- 
ing, it has greater flexibility in that it can be used in conjunction with 
both linear and nonlinear modulation and detection methods, with path 
selection on the basis of predetection as well as post detection signals. 
The principal diversity techniques would thus be space, frequency 
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or time diversity, in conjunction with "equal gain combining" or "selec- 
tion diversity." The error reduction afforded by the two latter methods 
is discussed below. 

7.8 Error Probabilities with Equal Gain Diversity 

The error reduction afforded by equal gain diversity transmission has 
been determined by Pierce 8 for binary FM with coherent and nonco- 
herent dual filter detection, on the premise of sufficiently slow flat 
Rayleigh fading, such that errors from noise alone need to be considered. 
For binary PM with differential phase detection, the error probability 
with equal gain diversity transmission has been determined by Voelcker, 
considering both errors from noise [P e (3) ] and errors from time variations 
in the transmittance [P e (2) ]. Voelcker has also determined the error 
probability with dual diversity transmission for four-phase modulation 
with differential phase detection, considering errors from transmittance 
variations with time alone. For all of the above cases, the following 
approximation applies for the probability of single digit errors with 
dual diversity transmission over independently fading paths 

P e , 2 ^3Pe,l (171) 

where P e ,i is the error probability for transmission over a single path 
(no diversity). For four-phase modulation, Voelcker's more exact 
expression, when reduced to small error probabilities, gives a factor 
4ir(3 + ir)/(2 + tt) 2 « 3.13 in place of 3 in (171). 

The mechanism responsible for error reduction by diversity trans- 
mission in the above cases also applies to transmission over channels 
with selective fading when the errors are caused principally by inter- 
symbol interference. With independently fading transmission paths 
there will be no correlation between intersymbol interference in the 
various channels, even though the signals are the same. Hence relation 
( 171 ) would also be expected to apply for the combined error probability 
P e given by (160). 

For small error probabilities, the following approximate expression 
is given by Pierce 8 for the error probability owing to noise with flat 
Rayleigh fading for binary FM and multidiversity transmission 

P„M ( ??~ 1 il. J, »' (172) 

ml (m — 1)1 

P,. 2 w3P.,i 2 (173) 

P e , 3 ft* 10/ 5 ,,! 3 (174) 

P B , 4 « 35P,,i 4 . (175) 
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The optimum number of diversity paths will depend on a variety of 
considerations, among them the available bandwidth and transmitter 
power, system complexity, and the source of errors. When the errors 
are caused by noise it is possible to realize a certain minimum total 
average signal power for a specified error probability P e , m , by appro- 
priate choice of m. As shown by Pierce and Harris, the minimum 
total average signal power is attained for any specified error probability 
when m is so chosen that in each diversity channel p ftJ 3, or about 5 
db, for binary FM with dual filter noncoherent detection. The number 
of diversity paths required to realize the minimum total average signal 
power is rather large, and the signal power reduction that can be realized 
with more than four paths is fairly small. For example, Pierce 18 shows 
that for an error probability P ei „, = 10~ , the minimum average signal 
power is realized with m = 16, for which the total signal-to-noise ratio is 
16.7 db, corresponding to a signal-to-noise ratio per channel of 4.7 db 
(p" = 2.95). With m = 1 the average signal-to-noise ratio is 40 db and 
with m = 4 is 19.4 db. Hence only a i-mall additional reduction in signal 
power is realized when the number of diversity paths is increased from 
m = 4 to m = 16. 

7.9 Error Probabilities with Selection Diversity 

Equal gain diversity as considered above entails a linear addition 
of the baseband outputs of the various demodulators, and would be less 
effective in conjunction with nonlinear demodulation methods, such 
as binary FM with frequency discriminator detection. With the latter 
method, switch or selection diversity reception would probably be pref- 
erable, in which only the receiver having the largest signal is selected. 
With this method the following relations apply for m-diversity transmis- 
sion when the errors are caused by noise and when receiver selection is 
based on the largest carrier signal at the detector input 

P e , m « 2 m - 1 m!P c .r (176) 

P c o ^4P C|1 2 (177) 

P e , 3 ^24P e>1 3 (178) 

P eA ^192P frl 4 . (179) 

For equal error probability, the average signal power with selection 
diversity must be greater than with optimum diversity by a factor 
equal to the mth root of the ratio of the factors in (176) and (172). The 
power must thus be increased by 0.62, 1.27 and 1.85 db for m = 2, 3 
and 4, respectively. 
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7.10 Multiband Digital Transmission 

The curves in Figs. 21 and 22 suggest that for a given total transmitter 
power and channel bandwidth, the error probability can be reduced by 
transmitting at a slower rate over each of a number of narrower channels 
in parallel. An approximate optimum bandwidth for each channel would 
be such that P e a) + P e (2) is minimized. This can be accomplished with 
separate transmitters and receivers for each channel, such that mutual 
interference between channels is avoided. Hence the adverse effects of 
selective fading can be overcome with the aid of more complicated 
terminal equipment, without the need for increased signal power or 
channel bandwidth. 

An alternative method that is simpler in implementation is to trans- 
mit the combined digital wave from the parallel channels by frequency 
or phase modulation of a common carrier, as ordinarily used for trans- 
mission of voice channels in frequency division multiplex. This method 
entails some mutual interference between channels, as well as greater 
channel bandwidth and carrier power than with direct digital carrier 
modulation, as discussed below. 

With the above method, the spectrum of the modulated carrier wave 
will have greater bandwidth than with direct digital carrier modulation. 
To avoid excessive transmission distortion of the combined wave, the 
bandwidth between transmitter and receiver must be at least twice 
that with digital carrier modulation. Hence, at least 3 db greater average 
carrier power is required in order that the noise threshold level of the 
common channel be comparable with that of direct digital carrier modu- 
lation. 

With such multiband transmission, intersymbol interference owing 
to selective fading is avoided, in exchange for mutual interference be- 
tween the various channels owing to intermodulation distortion caused 
by selective fading. Such intermodulation distortion is dealt with else- 
where (this issue, part 2) for a modulating wave with the properties of 
random noise, which is approximated with a large number of binary 
channels in frequency division multiplex. The results indicate that 
under this condition intermodulation distortion will cause less trans- 
mission impairment than does intersymbol interference in direct digital 
transmission. Hence multiband transmission by common carrier modula- 
tion permits a reduction in error probability in exchange for at least a 
twofold increase in bandwidth and carrier power. However, this reduc- 
tion in error probability may be less than can be realized with direct 
digital carrier modulation in conjunction with a twofold increase in 
bandwidth and signal power with dual diversity. 
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Error probabilities in binary multiband transmission by frequency 
modulation of a common carrier are dealt with by Barrow 21 on the premise 
of slow fiat fading over the combined band, so that only errors owing to 
noise need be considered and intermodulation distortion can be dis- 
regarded. 

VIII. SUMMARY 

The objective of this analysis has been to develop a transmission and 
modulation theory for troposcatter systems, applicable to digital trans- 
mission by AM, FM and PM at any speed and based on a realistic 
idealization of troposcatter transmittance properties. The basic model, 
together with the analytical procedure and certain basic assumptions, 
are reviewed here. 

8.1 Troposcatter Transmittance 

Based on certain physical considerations, an idealized multipath 
transmittance model is developed in which the received component 
waves vary at random in amplitude and phase and have transmission 
delays owing to path length differences which vary linearly with angular 
deviation from the mean path with maximum deviations ±A from the 
mean delay. With this type of model, a Rayleigh probability distribu- 
tion is obtained for the envelope of a received carrier wave in conform- 
ance with observations. 

To facilitate determination of transmission performance, two basic 
statistical parameters are required aside from the signal-to-noise ratio 
at the receiver. One of these is the autocorrelation function of envelope 
variations with time at a given frequency. The other is the autocorrela- 
tion function with respect to frequency at a fixed time. 

The first of these, the time autocorrelation function, depends on the 
rapidity of changes in the atmospheric structure within the common 
antenna volume. It has been determined by a number of observations 
with some theoretical support, as given in certain publications. 

The second basic parameter, the autocorrelation function with respect 
to frequency, has been determined by observation on a particular link. 
These observations conform well with the autocorrelation function 
determined analytically herein on the premise that the maximum delay 
deviation ±A noted above is given by the path length differences 
based on the beam angles between the 3-db loss points.* 

With the aid of this idealized model, endowed with the above basic 
parameters, as determined by observation or theory, it is possible in 

* This conclusion appears to be supported by the results of recent measure- 
ments on a 100-mile path. 2 * 
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principle to determine analytically the associated idealized transmission 
performance with any modulation method. Though an exact solution is 
possible in principle, it appears intractable and is not essential for 
engineering purposes. An approximate solution for transmission at 
any digital rate is derived herein. To this end certain basic statistical 
parameters are determined from the above two autocorrelation func- 
tions. 

8.2 Variations ill Transmittance with Time 

In Section II, distributions arc given for the time rate of change in 
the envelope and for the first and second derivatives of the phase func- 
tion. These probability distributions permit approximate evaluation of 
changes in the envelope, phase and frequency over a signal or pulse 
interval for narrow-band signal spectra. 

8.3 Variations in Transmittance with Frequency 

The corresponding probability distributions with respect to varia- 
tions in transmittance with frequency are given in Section III and permit 
approximate determination of random attenuation and phase distortion 
over the band of the signal spectra owing to the selectivity of fading. 
From these random variations it is possible to determine the correspond- 
ing pulse distortion together with resultant intersymbol interference in 
carrier pulse trains and error probability in the absence of noise. 

8.4 Errors from Selective Fading 

As a next step in the determination of error probability, an approxi- 
mate evaluation is made in Section IV of the probability of errors from 
intersymbol interference with selective Rayleigh fading in the absence 
of noise. In a first approximation it turns out that attenuation distortion 
can be neglected in comparison with phase distortion. Furthermore, the 
latter can be approximated by a component of quadratic phase distor- 
tion, or corresponding linear delay distortion. Intersymbol interference 
owing to quadratic phase distortion is determined for various carrier 
modulation methods, and an approximate relation is derived for the 
resultant error probability in the absence of noise. 

8.5 Errors from. Nonselective Rayleigh Fading 

With transmission at sufficiently slow rates, errors can occur in the 
absence of noise, owing to changes in amplitude, phase or frequency over 
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a pulse interval, caused by nonselective Rayleigh fading. The proba- 
bility of errors on this account is determined in Section V on the approxi- 
mate basis that changes over a pulse interval are proportional to the 
time derivatives of the amplitude, phase or frequency, depending on the 
modulation method. Comparison with available exact solutions for 
phase modulation shows that the inaccuracy resulting from this approxi- 
mation is inappreciable. 

8.6 Errors from Random Noise 

In Section VI expressions are given for the probability of errors from 
random noise with flat Rayleigh fading, as derived in various publica- 
tions for different digital carrier modulation methods. In addition, an 
expression is derived for error probability with rapid Rayleigh fading 
in conjunction with slow log-normal fading, as encountered on tropo- 
scatter links. 

8.7 Combined Error Probability 

In the final Section VII the combined error probability is determined 
on the approximate basis that it is the sum of the error probabilities for 
the three basic sources assumed above. Charts are presented from 
which can be determined the approximate combined error probabilities 
for binary phase and frequency modulation over a single path, and 
approximate expressions are given for the error probability with diversity 
transmission over independently fading paths. 

8.8 Basic Approximations 

The idealized model of troposcatter transmission assumed herein is of 
course an approximation, as are the idealizations regarding the per- 
formance of the carrier modulation methods. Even with exact mathe- 
matical analysis based on this model, the predicted performance would 
not conform entirely with that observed on actual systems. 

In determining error probability from the idealized model, two basic 
approximations were used to obtain numerical results. One is that the 
maximum departures ±A from the mean transmission delay can be 
determined from the beam angles taken between 3-db loss points. On 
short links with narrow-beam antennas, these are virtually equal to the 
free-space antenna beam angles, but for long links are greater owing to 
beam broadening by scatter. The second approximation is that errors 
from distortion owing to selective fading are caused principally by a 
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quadratic component of phase distortion. This is the first component 
that gives rise to distortion in a power series expansion of a nonlinear 
phase characteristic as a function of the frequency from the carrier. 

The same two basic approximations have been used in a companion 
paper (this issue, part 2) in a determination of intermodulation noise in 
analog transmission by FM of signals with the properties of random 
noise. Theoretical predictions based on free-space beam angles are in 
this case in reasonable agreement with measurements on two tropo- 
scatter links 185 and 194 miles in length, with narrow-beam antennas. 
Measurements on links 340 and 440 miles long give intermodulation 
noise that would correspond to beam angles and maximum delay differ- 
ences ±A that are greater than for free space by factors of about 1.35 
and 2.15, respectively. 

The above measurements also show that as the bandwidth increases, 
actual intermodulation noise will be progressively smaller than predicted 
on the premise of quadratic phase distortion. Translated to digital 
transmission, the error probabilities P« (1) owing to selective fading as 
determined here on the premise of quadratic phase distortion would 
represent an upper bound, that should conform well with actual error 
probabilities when the latter do not exceed about lO" 2 in Figs. 21 and 22. 

8.9 Comparison with Recent Related Publications 

Since the completion of the galley proof of this paper an article by 
Bello and Nelin 22 has appeared, dealing with errors in binary transmis- 
sion owing to frequency selective fading by a different analytical pro- 
cedure than used here. Numerical results are presented for error prob- 
abilities in dual and quadruple diversity transmission by binary FM 
with dual filter incoherent detection and binary PM with differential 
phase coherent detection. These results are based on an assumed Gaus- 
sian correlation function, or power spectrum, of the selectivity of fading 
with frequency. A comparison is made below of the above numerical 
results with those obtained on similar premises from relations presented 
here. 

For a Gaussian power spectrum of correlation bandwidth B c as used 
in the above paper, the corresponding value of a 1 in (18) is a 2 = 2(wB c )- 2 . 
Expression (55) applies with fo 2 /bo = a 2 in place of A 2 /3. With this sub- 
stitution and with T = B~\ expression (101) and Fig. 17 apply, with 
A- B = 0.79(B c T)-\ where {B c T)~ l is the parameter appearing in Figs. 5 
and 9 of the above paper for the irreducible error probabilities. 

Binary FM with dual filter detection as assumed in the above paper 
can be considered equivalent to ideal complementary binary AM over 
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each of two channels. When the frequency selectivity of fading is suffi- 
cient to cause errors in one or the other of these channels, the above 
method is essentially equivalent to dual diversity transmission by AM 
over two independently fading channels. On this basis, binary FM with 
dual diversity and dual filter noncoherent detection is approximately 
equivalent to binary AM with quadruple diversity. The error probabil- 
ities determined on the latter premise with AB = 0.79(B c T)~ l in (101), 
or in Fig. 17, in conjunction with (172) for m = 4, conform reasonably 
well with those given in Fig. 5 for dual diversity with \p = and n = 1. 
Complete agreement is not possible for the reason that the results in 
Fig. 5 assume a rectangular shape of undistorted pulses, whereas the 
present analysis is based on a more realistic pulse shape with a raised 
cosine spectrum, as indicated in Fig. 13. 

In the case of binary PM with differential phase detection, the rela- 
tions presented here with A -5 = Q.7§{B C T)~ X yield error probabilities 
that are significantly smaller than those given in Fig. 9 of the above 
paper. This is to be expected, since the present relations are based on 
detection with an optimum threshold level, whereas those in the above 
paper assume zero threshold, which is not the optimum owing to the 
presence of a substantial bias component in the demodulator output, 
when pulse distortion is pronounced. 13 Moreover, the shapes of the un- 
distorted pulses are different, as noted above. 

It is evident from the above considerations that apparently unrelated 
and possibly misleading results can be obtained unless comparisons are 
made of binary modulation methods of equal bandwidths with optimum 
implementation of each, as was done in Fig. 17. 

The above article called attention to another paper 23 by the same 
writers that refines Voelcker's original analysis 9 of errors in transmission 
over narrow-band channels owing to transmittance variations with time. 
Their results show that for a Gaussian power spectrum of the fading 
rate as assumed herein, Voelcker's analysis is exact, though this is not 
true for all forms of power spectra. 
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APPENDIX 



Transmittance of Troposcatter Channels 

Owing to the differences in path length from transmitter to receiver 
via the various heterogeneities in the common volume, the various 
components of the received wave arrive with different delays. For 
analytical purposes it is convenient to assume a certain mean reference 
path with delay T and to express the transmission delay via other 
paths relative to the delay T . Actually there will be a large number 
of paths with the same delay T as the mean path and a large number 
of paths for each other delay. In the present analysis the approximate 
model indicated below is assumed, with a single vertical scatter plane 
midway between transmitter and receiver. 




The amplitude of the wave component arriving over a path at the 
distance x above the mean path is taken as A(x,t) and the delay over 
this path as 

T(x) = To + S(x). 

The wave component arriving via this path is then 

e x (o},t) = A(x,t) cos u[t - To - «(*)]. (180) 

Let L be the distance between transmitter and receiver and H the 
height of the mean path. In this case 



8(x) = s(x)/v 



(181) 
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where v is the velocity of propagation and s(x) the path length difference 
given by 



t(x) =[|+(H + a;) 2 J-(f + H 2 ). 



(182) 



(185) 



In actual systems H « L. Furthermore, the maximum value x of x is 
ordinarily much smaller than H. On these premises the following ap- 
proximation applies 

50) = (2H/Lv)x = x/c (183) 

where c = vL/2H. 

It will further be assumed that there is an infinite number of paths, 
in which case the received wave becomes 

e(w,0 - / A(x,t) cos u(t - To - x/c) dx (184) 

= cos u(t - To) \ [A(x,t) + A(-x,l)] cos (wx/c) dx 

4- sinco(i - To) /"* [A(x,t) -A(-xfi)]mn (ax/e) dx. 

It will now be assumed that 

f X [A(x,t) + A(-x,t)]dx = 0. (186) 

■'o 

This appears to be an appropriate physical requirement, for the reason 
that reflections occur as a result of variations in the electrical properties 
of an elementary volume, relative to that of the common volume. No 
reflections occur with a uniform common volume. In a heterogeneous 
common volume, each positive reflection must be accompanied by an 
equal negative reflection, which is reflected in condition (186). More- 
over, under this condition there is no reflection along the mean path of 
the transmitted beam. That is, with x = in (185), e(t) = provided 
(186) applies. 

Condition (186) can be insured if the following Fourier series repre- 
sentations are used for x ^ x 

A{x,t) + A(-x,t) = J^a(m,t) cos niTx/A (187) 

m=l 

and 

A(x,t) - A(-x,t) = E&(w,0 sin mwx/x. (188) 
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With m = 1, 2, 3, etc., as above, the area under each harmonic com- 
ponent vanishes, such that condition (186) is satisfied. 

With (187) and (188) in (185), the following relation is obtained 

c(w,0 = cos u(t - T)U(<a,t) + sin <a(t - T)V(<a,t) (189) 

where 

U(u,t) = 2 a(m,t) / cos rrnrx/x cos cox/c dx (190) 

m-l •'0 

V{ta t t) = 2-i b(m,t) / sin rrnrx/x sin cox/c dx (191) 

Evaluation of the integrals yields the following expressions 
U(«,t) = ± \A( m ,l) f »i° (mr - «A) Bin (m, + ^n () 

m=l L «T — CoA W7T + CoA J 

F(»,0 = ± IBM [~ sin (mr " » A) - sin (TOr , ± " A) 1 (193) 

m=l L m ^ ~ ^A W7T + CoA J 

where 

A(m,t) = £a(m,t) 

B(m,t) = &b(m,t) (194) 

A = x/c. 

It will be noted that A is the maximum departure from the mean 
delay T . 

In evaluation of (192) and (193) it is convenient to introduce a new 
reference frequency coo in place of 0, and to choose this reference fre- 
quency such that 

cooA = n-K. (195) 

Thus 

uA = tit + wA (196) 

where — t < wA < w, and u is the deviation in frequency from coo . 
The functions (192) and (193) are then replaced by 

U(u,t) = 2^ \Mm,t) <— -^ ^ — - J 

m=i I (m — n)T — uA 

(197) 
, sin [(w + n)ir -f- wA]l 

sin (m + u)t + wA 
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tri ,\ V ib/ ,\ /sin [(w — n)ir - i<Al 
STi i (m — n)w — uA 

_ sin [(ro + n)w + mA] \ 
sin (m + ?i)tt + mA J 



:i98) 



Iii troposcattcr transmission it turns out that m is of the order of 
100 to 1000. For this reason the second terms in the above series, in 
(m + n)ir, can be neglected. With this simplification and with m — 
n = j, expressions (5) and (fi) are obtained. 

Expression (189) can then be written in the form 

e(u,t) = r(u,l) cos [u(t - T) - <p(u,t)] (199) 

where r and v are given by (3) and (4). 

The channel transmittance is accordingly given by (2). 
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